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corrected vertices, which are applicable in general, and provide a recipe for the inclusion 
^ ' of the dominant and subdominant BLO effects in existing LO calculations. We also derive 

' similar expressions in the mass insertion approximation. We investigate in more detail 

the focusing effect pointed out in our earlier work, which, at large tan/3 and /U > 0, leads 
to a reduced supersymmetric contribution to the above processes. We also find that, in 
some cases, flavour dependence, that accidentally cancels at leading order, can reappear at 
BLO. We further include electroweak corrections, which, while generally subdominant, in 
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1. Introduction 

Flavour physics, in both the leptonic and hadronic sectors, currently provides one of the 
best hopes of discovering, or at least constraining, new physics beyond the Standard Model 
(SM). In the hadronic sector in particular, decays mediated by flavour changing neutral 
currents (FCNC) play an important role as the Glashow-Iliopoulos-Maiani (GIM) mech- 
anism Q ensures that both SM and contributions due to beyond the SM (BSM) physics 
enter at the one-loop level. It is therefore possible that such contributions can be com- 
parable to the SM ones, or even completely dominate the behaviour of the underlying 
process. Once one takes into account the increasingly accurate experimental data that is 
being gathered at both dedicated flavour physics experiments, as well as the B~physics 
programmes operating at collider experiments, useful constraints can often be placed on 
the parameters and mass scale of a given model of new physics. Conversely, for some 
processes, like Bg — > ix'^ ^~ , a measurement of a branching ratio at the Tevatron would 
immediately indicate a detection of BSM physics. 

One of the most compelling extensions of the Standard Model is the Minimal Supersym- 
metric Standard Model (MSSM) |^. The non-renormalization theorem of the underlying 
supersymmetric theory can explain the stability of scalar potentials in theories involv- 
ing two different hierarchies. Additionally, the MSSM provides a viable cold dark matter 
candidate (namely the lightest supersymmetric particle), a natural scheme for gauge cou- 
pling unification, and is usually compatible with the precision electroweak data currently 
available. 

Softly broken low-energy supersymmetry (SUSY), however, like most new physics 
schemes, allows for the possibility that contributions to FCNC and CP violating processes 
can exceed SM expectations by orders of magnitude (the flavour and CP problems). The 
source of the flavour problem in the MSSM is primarily due to the arbitrary nature soft 
supersymmetry breaking terms ||2|. 
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The most common approach to these problems is to assume that the underlying the- 
ory obeys the conditions imposed by minimal flavour violation (MFV) The definition 
of MFV, presented in is that flavour violation is determined completely by the struc- 
ture of the usual Yukawa couplings. In other words, the mixings among the down and 
up squarks are governed by the CKM matrix. In the MSSM this restricts the form the 
soft terms can take (see for the exact expressions). One popular scheme, that respects 
MFV, is that the soft terms are universal at some high scale associated with supersymme- 
try (SUSY) breaking, like the grand unified scale or the Planck scale (a parameterisation 
used, for example, in the Constrained MSSM). However, this hypothesis is not renormal- 
ization group invariant. Flavour violating terms are induced, via running from the high 
scale A to the characteristic mass scale of the squarks ususY, that are proportional to 
log(A^///|[j^y)/(47r^) IQ]. It should be noted however, that, provided the theory satisfies 
MFV up to the scale A (a rather strong assumption), all FCNC transitions remain propor- 
tional to the appropriate CKM matrix elements and the resulting low energy theory still 
satisfies the MFV hypothesis. However, once seeds of non-universality are introduced at 
the high scale it is possible that they can become amplified by running. 

This provides motivation to generalise to a broader framework, namely general flavour 
mixing (GFM) in the sfermion sector. In general, the flavour structure of the soft terms is 
not protected by any symmetry and can be rather arbitrary. One simple example is that 
a degree of non-universality can be allowed in the squark soft terms (beyond that allowed 
by MFV). In this case additional effects are possible that are proportional, essentially, to 
the degree of splitting between the entries for each generation. 

Deviations from MFV can easily appear in a variety of SUSY models. In theories 
with SUSY breaking mediated by supergravity, for example, it is possible to induce a 
wide range of flavour violating effects Q once one proceeds beyond the simplest minimal 
SUGRA scheme Grand unifled theories involving right handed neutrinos, like the 
minimal SO(IO) models with a speciflc family structure, often lead to additional sources 
of flavour violation due to the interactions that exist, at the unification scale, between 
right-handed down squarks and neutrinos ^, ^, ^ |l^] . Experimental limits and results 
are therefore especially helpful when restricting the possible mixings between the various 
generations and constraining these models. 

In this paper we shall concern ourselves chiefly with flavour violation between the sec- 
ond and third generations. FCNC processes involving such transitions have been studied 
in detail in the context of the SM and the short-distance contributions to a wide variety of 
processes have typically been calculated to NLO in the SM (the evaluation of long-distance 
effects is another matter). In the case of B ^ Xs'j these efforts have resulted in a very 
good agreement between SM calculations and experimental results with relatively little 
room left for new physics. When placing constraints on a given model it is useful to have 
a calculation that is of a similar accuracy to the SM contribution. In the MSSM com- 
plete NLO calculations, however, are rather complicated as additional two-loop diagrams 
involving gluinos need to be evaluated. It is, however, possible to include the effects that 
are large once one proceeds beyond the LO (BLO). Such effects are typically classifled as 
being proportional to either tan/5 or large logs. Such BLO analyses have been performed 
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in MFV ll, H, |, ll, |1J and, more recently, in GFM |l5|, y]. In GFM, in partic- 
ular, a focusing effect was found in [l^ that gave rise to significant shifts in the allowed 
regions of parameter space compared to a LO analysis. (A similar effect appears also in the 
MFV scheme but is much weaker |l^.) Basically, in many cases of phenomenological 
interest {e.g., large tan/3 and > 0), SUSY contributions to B ^ Xg^ are significantly 
reduced compared to the LO approximation. A similar effect was also found in the decay 
Bs fJ'^fJ'~ and Bs — Bg mixing [p!7|]. 

The aim of this paper is to present the first complete analysis of dominant BLO effects 
in general fiavour mixing in the case of three processes. The first, B ^s7, has been 



discussed previously in |15, 16 1. However, here we shall include the additional corrections 
that arise when one includes charginos and neutralinos in the resummation procedure 
discussed in [|^]. In particular, we include contributions arising from higgsino exchange, 
that are proportional to the Yukawa couplings of the third generation, and the additional 



contributions to the charged Higgs vertex that were discussed in the context of MFV in |14|. 

The other two processes we shall consider are the decay Bg /U^/t^^ and Bg — Bg 
mixing. These processes have not been observed yet but have come under a lot of theoretical 
scrutiny lately due to the large contributions possible in the large tan/3 regime. In this 
paper we discuss the GFM contributions to both processes in detail, highlighting the effects 
that appear once one proceeds beyond the LO. 

In all the three cases, we shall present the analytic expressions required to implement 
BLO corrections in the GFM scenario for possibly large deviations from the MFV scheme. 
However, since these general expressions are often rather complicated, we shall also derive 
expressions in the the mass insertion approximation (MIA), allowing the BLO effects to 
be shown explicitly. In both cases we will provide an explicit recipe for including the 
BLO effects into the existing LO expressions. Whilst we shall not include such effects 
in the forthcoming analysis, the formalism we shall present should, with relatively little 
modification, be applicable to the GP violating case. 

The paper is organised as follows. In section || we summarise the formalism employed 
in this paper, giving complete expressions for all the corrected masses and vertices used in 
our calculation. In section |^ we present analytic expressions for these masses and vertices 
in the MIA. In section |^ we discuss the decay B — > Xg'j providing analytic expressions 
for the BLO corrections to supersymmetric and electroweak contributions in the MIA. In 
sections ^ and ^ we perform a similar analyses for the decay Bg — > fi^ and Bg — Bg 
mixing, respectively. Finally, in section ^ we present our numerical analysis. 

2. Beyond Leading Order Effects and General Flavour Mixing 

The influence of tan /3 enhanced effects on the down quark masses, the charged Higgs and 
neutral Higgs vertex are known to be large. It is therefore essential, especially when working 
in the large tan (3 regime, that such contributions are taken into account (and resummed 
if necessary). 



In this section, we shall follow the method flrst developed in |15, and generalise it to 



include the additional effects that appear once the contributions of chargino and neutralino 
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loops are taken into account. It should be noted that the analysis below encompasses both 
MFV and the GFM scenario and can be easily extended to include, for example, CP 
violation or flavour violation in the leptonic sector. 

2.1 The Framework 

Once the supersymmetric particles have been integrated out, the effective Lagrangian de- 
scribing the quark mass terms, at some scale ^ < MsjjSY, in the physical super-CKM 
basis (SCKM) is given by 

-C^^' = dR (mf + 6md) dL + UR (mi") + Sm^) ul + h.c, (2.1) 

where dL,R and ul,r denote the down and up components of the left and right quark fields, 
respectively.^ In the physical SCKM basis the quark mass matrices are, by definition, 
diagonal and it is possible to make the identifications 

md ='rn~d^ + 5md = diag (m^, ms,mb) , (2.2) 
niu =m^'> + 6mu = diag {rriu, rric, mt) , (2.3) 



where md,s,b and mu,c,t denote the physical masses of the down and up-type quarks re- 
specitively. The bare mass matrix m^^^ is rela 
from the superpotential in the usual manner 

i,u ~ d,u 



specitively. The bare mass matrix rn~^'^ is related to the 3x3 Yukawa couplings Y^^^ derived 



4°i = -^,«4°i- (2-4) 



where Vd^u = i^dv)- Finally, 6md and 6mu denote the radiative corrections to the quark 
masses induced by integrating out the SUSY particles [18, 19, 20, ^]. The corrections have 
the form^ 

6md = Y^L + )fiRmf + (2.5) 

(5m„ is given by a similar formula after one performs the substitution d ^ u. The 3x3 
hermitian matrices S^']^ ^ and the 3x3 complex matrices S^*^ denote the contributions 
arising from wavefunction and mass corrections due to two point diagrams involving gluinos, 
charginos, neutralinos and squarks. (Full expressions will be given later in the text.) 

Before discussing how the radiative corrections brrid are calculated, it will be useful to 
consider the transformation from the interaction basis to the physical SCKM basis. In the 
interaction basis, the MSSM superpotential is 

Wf = -tiHdHu + Yl^^^HdL'^E'' + YP°HdQ''D° - Y^'^^HuQ^U", (2.6) 

Q" and L° are the quark and lepton SU(2) doublet superfields, D°, U° and denote the 
singlet superfields and Hu and Hd are the two Higgs doublets that appear in the MSSM 



^It should be noted that, throughout this section, we shall adopt matrix notation and suppress flavour 
indices unless otherwise specified. 

^As we allow the inclusion of electroweak effects we will not assume proportionality to the strong coupling 
constant here, unlike in llsL Iql. 
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(for more details see, for example, [||), while Y°^^ are the appropriate 3x3 Yukawa mass 
matrices in that basis. 

The soft SUSY breaking terms are also usually introduced in the interaction basis. 
The Lagrangian for the bilinear soft SUSY breaking terms is given by 



/^o mass 

'"q soft 



+ i'^mldl + 4^m2,4 + [d"^ {vdA*d) 4 + h.c. 

I -of 2 ~o I -of 2 ~o I ~o\ t A*\ ~o I u 



(2.7) 



where mg, and mfj are, in general, arbitrary 3x3 hermitian matrices. The trilinear 
terms and A^, on the other hand, are arbitrary 3x3 complex matrices. We have not 
assumed that the trilinear soft terms are proportional to the appropriate Yukawa coupling. 
(We discuss an alternative parameterisation, that can be used in the GFM scenario, in 
appendix p.) 

Transforming the quark fields from the interaction basis to the physical SCKM basis 
involves performing unitary transformations on both the left and right handed fields such 
that 



dR = ydRd°R 



dh = Vd^dl, 



UL = Vii.r dr. 



(2.8) 
(2.9) 



The bare mass matrix is then related to the Yukawa couplings defined in the interaction 
basis via the relation 



m 



(0) 



™{0) _Y y(0)oyt 



(2.10) 



m 



(0) 



and rriu^ appear in all quantities derived from the superpotential (2^) not subject to 



the corrections ( |2.5D , such as the couplings of supersymmetric particles. The CKM matrix 
K is related to the transformations ( |2.8| )-( ^^ in the usual manner 



(2.11) 



As the radiative corrections SrridM are calculated in the SCKM basis, it is necessary to 
consider how the transformations ( |2.8D -(^l9|), when performed on the squark fields, affect 
the relevant mass matrices. After transforming to the physical SCKM basis, the soft terms 
become 



m 



III =^i3L™Q^Di' 

2 T/t 



,RR =VDRmj)Vl^, ml^j^ =Vdi^ {vdA}) V^^, 



•mtLL =yUi^rin?QV) 



ml^RR =VuRmuV^^, 



mlLR=VuAvuA:)V^. (2.12) 



The 6x6 down squark mass matrix may then be written in the following manner 



d 



^IlL + Pd,LL + Dd,LL ^2 + Fd,LR 



^d,LR + ^d,LR 



^d,RR + ^d,RR + Dd^RR 



(2.13) 
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(The up squark mass matrix may be similarly defined by substituting d with u.) The 
F-terms that appear in (2.13) are given by 

i^d,LL = t^i°\ i^d.iji? = t, Fd,LR = -/xtan/3mi°)t, (2.14) 

and the flavour diagonal D-terms are 

Dd,LL = mz cos 2/3 {T^d - Qd sin^ Ow) h, Dd,RR = mz cos 2[3Qd sir? Ow'^'i- (2.15) 

It should be noted that, in the physical SCKM basis, the F-terms are, in general, not 
necessarily flavour diagonal as they are derived from the superpotential and are therefore 
functions of the bare mass matrix ma^P . 

To obtain the physical squark masses m^^ it is necessary to perform an additional 
unitary transformation on the squark fields such that 

<2rt 



r,A^4rT=diag(^m,-^,...,m,-J. (2.16) 

It is conventional to decompose the original 6x6 unitary matrix Yd into two 6x3 subma- 
trices T^l and TdR: 

{Td)u={TdL)u^ (rd),,+3 =(rrfH)/.. (2.17) 

where I = 1, . . . , 6 and i = 1, 2, 3. 

Departures from the MFV scenario are often parameterised in terms of the dimension- 
less quantities 



''^d,RL I .. / , \ ( ^d,RR , 

Srl)^.=^^==^=, {4r)^.=^= ' (2.19) 

' ^^Irr) . . (^Ill) . . J (^Irr) . . (^Irr) . . 

The soft terms j^y Y = L,R) are given in ( |2.12| ) and i,j = 1,2,3. Similar definitions 
apply for the up quarks. It should be noted that, since and j^j^ are related to one 

another by SU(2) invariance, we have the relation 

4l = K^IlK (2.20) 

Let us briefiy comment on the basis dependence of these definitions of Physical 
quantities such as cross-sections and branching ratios are naturally independent of the 
basis in which one defines the soft terms. The basis in which one defines the insertions 
(5^y, however, is essentially a matter of convenience. As discussed above, we work in the 
physical SCKM basis throughout this analysis and as such the definition ( ^.18 )-( ^.19 ) is 
essentially the easiest to implement numerically. Other definitions of d'^y have been used 
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in the literature before, for example, one might define S^y the uncorrected (bare) SCKM 
basis where the Yukawa matrices derived from the superpotential are diagonal (we shall 



define this basis more formally in subsection 3.4). Transforming between different bases 
involves performing additional unitary transformations on the soft terms ( |2.12| ) and, unless 
large non-universalities exist, the differences between the transformed and the original 
typically rather small. Below we will derive many expressions in the MIA where 
one usually assumes that the diagonal entries of the soft terms are universal. In light of 
the above we expect them to be applicable to alternative definitions of Sj^y- During our 
numerical analysis we employ a similar definition for the SUSY soft terms to ensure that 
our formalism remains applicable to as wide a variety of models as possible. 

After defining our framework, let us now move on to the effects that these corrections 
have on the electroweak and supersymmetric vertices. 

2.2 Corrections to Electrovi^eak Vertices 

Integrating out the supersymmetric particles, coupled with the effect of transforming be- 
tween the interaction and the physical SCKM bases, can affect the form of the electroweak 
(i.e. the Higgs and gauge boson) vertices present in the resulting effective theory. 

After transforming to the physical SCKM basis, the W boson vertex has the following 
form 

Cw = ULl'^C^dLW^ + URrC^dRW^ + h.c. (2.21) 
The 3x3 coupling matrices and , are given by 

= + IKlK + l^^iL^ + ACf = -^K^^ (2.22) 

= AC^. (2.23) 

We employ the notation AC^ to denote the vertex corrections that arise from three point 
diagrams when one integrates out the SUSY particles. Identifying the left handed coupling 
of the W boson with the physical CKM matrix K'^^^ that is measured from experiment, 
we have the relation 

K^^ = K + l^^^^K + iK^t^ - y^ACY- (2.24) 
2 2 c/2 



The uncorrected CKM matrix K is defined in ( ^.11 ) and appears in all vertices not subject 
to the corrections ( |2.24 ). 



Now consider the coupling of the Z boson with down quarks 

Czo = dLrcldLZl + dR^C^dnZ''^. (2.25) 
The 3x3 coupling matrices C£ and are given by 
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v2mvFsmp V ^ ^ 



The radiative corrections to C£ and can induce off-diagonal elements to the coupling 
that lead to additional sources of FCNC. 

Turning to the Higgs sector, the inclusion of radiative corrections is especially impor- 
tant. As the coupling between the Higgs sector and squarks features a dependence on the 
soft SUSY breaking terms (rather than only gauge interactions), the corrected vertices that 
arise when one integrates out the coloured SUSY particles can display a non-decoupling 
property. Large corrections to the vertices are therefore feasible for even TeV scale sparticle 
masses. 

Once one has integrated out the SUSY particles, the charged Higgs interaction becomes 
Cs+ = URCtdLS+ + ULCf^^dRS^ + /i.e., (2.28) 
where S'^ = , and the 3x3 matrices coupling are given by 

(2.29) 
(2.30) 

where = cos /3, sin /? and = sin /3, — cos /?. 

The neutral Higgs and Goldstone boson interact with the down quarks in the following 

way 

£50 = 4cf di^o + dLCtdRS\ (2.31) 

where = , hP , and the effective vertices and may be written in terms 
of the 3x3 matrices 

(2.32) 

and = cos a, — sin a, i sin /?, —i cos /?. 

In the limit where the physical SCKM basis is identical to the bare SCKM basis 
(i.e. where = rn'^'^), ( 2.32| ) is identical to the diagrammatic result derived in the on- 
shell formalism used in |23[ |. 

2.3 Corrections to Supersymmetric Vertices 

As the corrections to the electroweak vertices are calculated in the physical SCKM basis, it 
is necessary to discuss how the supersymmetric interactions are altered once transformed 
into this basis. Ignoring the effects of wavefunction renormalizations,^ that are not en- 
hanced by tan /3, the changes introduced by transforming to the physical SCKM basis 



^We do, however, include these contributions in our numerical analysis. 
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typically amount to the introduction of rn\^ and K into the various vertices. For instance, 



after these replacements have been performed, the chargino vertex becomes^ 

ft 

a.i.J 



= E 4 {x), [{CciL),j, Pl + {CdR)aM Pr] iA , (2.33) 



where o = 1, 2, i = 1, 2, 3, and J = 1, . . . , 6 and the couplings CdL and CdR may be written 
in terms of the matrices 

{CdL)aJ^ = - 92V:^ (F^L^) J, + ^ . ^ ^2 (r^ism^x) , (2.34) 



iCciR)aM -;U:2 [ruLKm^ , (2.35) 

\/2mw COS p ^ ' 



where K is defined in terms of the physical CKM matrix by K'^'^ in ( 2.24| ) and the bare 



masses are given by (^^) and its analogue for the up quarks. The matrices U and V 
diagonalise the chargino mass matrix ^A.^± such that 

\JM-^±V'^ = diag ("^^ij^^^i 

The appearance of the bare quark mass matrix Wb'^ in these vertices can lead to large 
effects in both MFV and GFM models. A full list of vertices relevant to our calculation 
appear in Appendix |C|. 

2.4 Numerical Aspects 

Let us now discuss how the method discussed above should be implemented numerically. 
As we will be investigating values of up to 0(1) for the flavour violating parameters (5^y 
(X, y = L,i?) (|2.18| )~( p.l9| ), it is important to devise a method such that the effects 
discussed in section ^ are taken into account, whilst also retaining the numerical accuracy 
associated with working in the squark mass basis. Such an iterative method was proposed 
in 1 16] and it will be useful for our purposes to briefly summarise it here. 



Once the unitary transformations ( |2.16| ) have been performed on the squark fields the 
gluino contribution to bm^ becomes |16| | 



^"^d = ^ ^2 (3) ^ {X*dR)u {^dL)ij rUg Bo 
7=1 



The Passarino-Veltman function Bq can be found in appendix |A.5| . Using this relation it 



is possible to calculate the bare mass matrix using (|2.3| ). It should be noted, however, that 
6md contains a dependence on m^^^ as it appears in the squark mass matrix through the F- 



terms ( p. 14 ). It is therefore necessary to employ an iterative procedure such that m^^^ and 



the mixing matrices TdL,R are determined to the desired level of accuracy. The inclusion 
of the effects induced by chargino and neutralino contributions introduces a dependence 



*Our notation for the supersymmetric vertices differs slightly from that used in |l^], broadly speaking, 
one may convert between the two by making the substitution L ^ R. 
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on and the uncorrected CKM matrix K in the formula for Srud- One must therefore 



expand and generahse the iterative procedure presented in such that these effects are 
included as well. 

In the first step of the procedure, iri~^^ , rriu^ and K are set equal to the input parameters 
rrirf, niu and K^^ , respectively, and 6md,u is set equal to zero. In the second step, the squark 
mass matrices and the super symmetric couplings are then calculated with these input 
values, allowing the evaluation of the radiative corrections S"^, T,^ and ACff^ using the 
formulae presented in appendix |D|. In the third step, Snid^u: the bare mass matrices and K 
are determined using (2^) and ( ^.24 ). The resulting values are then used as the new input 



parameters in step two. The second and third steps are then repeated until convergence 
occurs. The iterative procedure converges rather rapidly and after n iterations amounts 
to including the first n terms that arise in a Taylor expansion in tan /?. For an example 



of the procedure applied to MFV, see |16|. With the final forms of the supersymmetric 
couplings, uncorrected CKM matrix and bare mass matrices determined, the corrections 
to the Z boson, charged Higgs and neutral Higgs vertices may be calculated using the 



formulae presented in section 2.2 and appendix 0. We should emphasise here that we work 



in the squark mass basis throughout and therefore include all the effects that can occur at 
higher orders in the MIA as well as the BLO effects described in the previous subsections. 
In addition, we include the effects induced by additional electroweak contributions, light 
quark effects and SU(2)i x U(l)y breaking. 

3. The Mass Insertion Approximation 

,{o) 



Expressions for and the corrected vertices are well known in MFV models |2^, 11, 
and it will be useful for our purposes to extend these results to the GFM scenario. The 
flavour dependence of analytic expressions is often rather obscure when flavour violation 
is communicated via the matrices ( |2.17| ). To express the underlying dependence on the 
off-diagonal elements of the soft breaking terms it is therefore useful to work in the mass 
insertion approximation (MIA) . According to this approximation the off-diagonal elements 
of A^j are treated as perturbations and enter expressions through mixed propagators pro- 
portional to the relevant element (or insertion). These insertions are parameterised in terms 
of the dimensionless quantities defined in ( |2.18| )-( 2.19| ). Equivalently, one may expand the 



matrices ( |2.17| ) about the diagonal. When performing actual numerical calculations it is 
more advantageous to diagonalise the squark mass matrices ( |2.13| ) using numerical routines 
to ensure that higher order terms in the MIA are included, this is what we shall do in our 
numerical analysis presented in section |^. 

Before proceeding with our analytic expressions for the bare mass matrix and the 
various effective vertices, let us first outline the approximations we shall use throughout 
this section. As we are chiefly concerned in exhibiting the dominant behaviour displayed by 
the corrections to the bare mass matrix and effective vertices, in this section we shall work in 
the approximation of vanishing electroweak couplings and typically ignore SU(2)i x U(l)y 
breaking effects. We therefore mainly concern ourselves with the effects induced by gluino 
exchange, and those that arise from higgsino exchange that are proportional to the Yukawa 
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/ di^ dii \ I til ii \ 
di. ' ^ dp, di. < ^ dr. 
— ■ — mT6m61515m^ - — — • — ^ — 

Figure 1: The dominant gluino and higgsino diagrams tliat contribute to S^^. 



couplings of the third generation. Let us emphasise, however, that during our numerical 
analysis we include all the effects that arise from non-zero electroweak couplings, S\]{2)l x 
U(l)y breaking effects and the effects of the Yukawa couplings of the first two generations. 
Concerning the accuracy that we work to within the MIA, we typically include terms up 
to second order in the MIA (unless specified otherwise). We therefore do not include the 
effects of multiple diagonal LR insertions that are proportional to lr~ f^''^d!^ tan/3. It is 
possible to resum the effects induced by these insertions via the method outlined in |23]. In 
this method however, some BLO effects are then encoded into the factors of cos 6^ and sin ^ 
that appear in the MFV squark mixing matrices. Since the ultimate aim of this section 
is to present analytic expressions that represent all of the BLO corrections that appear 
in the framework presented in section |2|, we shall only consider the effects of at most one 
flavour diagonal LR insertion. Converting our expressions to take into account such effects 
however should be relatively easy. Finally, to allow for easy comparison between our results 
and those that already exist in the literature, we have assumed that the trilinear up-type 
soft terms are proportional to the appropriate Yukawa coupling (i.e. = m^Au, where 

Au = diag{Au, Ac, At)), although we still do not assume a similar relation for the down- 
squark sector, (it is easy to convert back by making the substitution At — > (jn^ ijij /nit 
in the various expressions that follow). 



3.1 The Bare Mass Matrix in the MIA 



We shall first consider the corrections induced by CKM and GFM effects on the bare mass 
matrix m^^\ Loop corrections to the bare mass matrix in the MFV scenario were derived 
in |ll|. They were subsequently generalised to GFM in |2^, 15, [T^ ]. 

The bare mass matrix may be determined by evaluating the self energy corrections 
that appear ( [2.51 ). The dominant contributions to S^^^ arise from self-energy diagrams 
involving gluino and chargino exchange, which are depicted in Fig. |l]. On the other hand, 
the corrections to ^ ^ are rather small when compared to those that arise from ^ 
as they are not enhanced by tan/3, nor do they feature a chirality flip on the gluino line. 
Coupled with the suppression factors of m^^^ that accompany them in (2^), their omission 
will not dramatically affect our final results. 



- 12 - 



Figure 2: Additional GFM contributions to the diagonal elements of 



(0) 



To first order in the MIA, the diagonal elements of m^^^ are given by 



m 



(0) 



1 + €{ tan j3 



(3.1) 



where i = 1,2,3 and denotes the combined dominant gluino and chargino contributions 



ei = es + evY^Si'i, 



(3.2) 



Yt is the top quark Yukawa coupling and is the usual Kronecker delta function. The 
coefficients and ey are given by 



ey = - 



At 



nin 



H2 {yuR^Vui) ■ 



(3.3) 
(3.4) 



In the above expressions is the strong coupling constant and 6*2(3) = 4/3 is the 
quadratic Casimir operator for SU(3) and At = {Au)^^- The loop function H2 is given 



in appendix A.l, whilst the arguments of the function are 



m 



Xj 
dL 



d,LL 



m: 



u,LL 



2 ' 



JUL 



(3.5) 



the definitions for x^^ and y^^^ may be easily obtained by substituting L with R in the 
above expressions. It should be noted that the soft terms that appear in ( ^.5| ) are common 
values of the diagonal entries of the SUSY soft terms ( p. 12 ). 

It is possible in GFM models to induce large contributions to the bare down and strange 
quark masses through diagrams involving three insertions (Fig. ||) |^6|. For example. 



m 



(0) 



'22 (1 + es tan /?) 
where €4 is given by 

64 = 



€4 tan j3x^^x^^ 



nis (l + eatan/?) 



'^^^32 V'^^«/23 



(3.6) 



(3.7) 



^In the following we shall neglect the flavour diagonal contributions that arise from the soft terms nn?d.LR 
unless they tan (3 enhanced. The corrections induced by these terms are included in our numerical analysis. 
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Figure 3: GFM contributions to the off-diagonal elements of 'm'"^\ 



The loop function H4 is given in appendix A.l whilst its arguments are given in ( |3.5D . 

Now let us turn to the off-diagonal elements of m^^^ . The diagrams in Fig. |^ and Fig. ^ 
illustrate the flavour violating corrections that arise from MFV and GFM contributions. 
Evaluating all four diagrams, we find the contribution 



V 1 + ej tan /? V ^^J ^J (1 + g . tan p f 

[eLLX^^ [5ii).. + eYY^Kt,Kt,] (m,),,tan/? 
(1 + eitan/3) (1 + ejtan/?) 

where erl and ^ll are given by 



(3. 



^RL = -^C2{3)H2 , CLL = -^C'2(3)^//3 [x^^, X^^, X;^^) (3.9) 

and Err can be obtained by making the substitution L <-> i? in the formula for cll- The 



loop functions H2 and are defined in appendix A.l, the CKM matrix K is defined 
m pUD , ei and ej are defined in ( |3.2| ), and x-j^^^ is given by 



^d,LL''^d,RR 



-s,. = ^-^^ • (3.10) 

It will be useful to see how the above expressions behave in the limit of degenerate 
sparticle masses. For instance, the various e-factors that appear in the above formulae 
become 

es=^s9n(^), =^^'^5?^ ( — ) , (3.11) 

37r 327r^ V / 



^G=^^sgn{\i) , eRL=^, eiL = - ^sgn {^) . (3.12) 



From ( p. 11 ) it is easy to see that, in the phenomenologically favoured region fj, > 0, Af < 0, 



the chargino and gluino contributions in ( [j.2| ) for i = 3 partially cancel. This can lead to 
a reduction of BLO contributions compared to a case where only gluino contributions are 
taken into account. 

As we are chiefly concerned with flavour violation in the down squark sector, we can 
safely omit the effects induced by LR, RL and RR mixings amongst the up squarks. In 
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Figure 4: The gluino correction to the Z boson vertex that arises at second order in the MIA. 



other words we assume that and are diagonal matrices. However, the insertion 

5^1^ is related by SU(2) symmetry to Sf^j^ and its effects on the bare mass matrix should 
be included. In the approximation used in this subsection however, the contributions 
proportional to S'^^, that arise solely from higgsino exchange, are rather small, as they are 
suppressed by factors of the Yukawa couplings of the first two generations. We will see in 



section 3.3 however, that, once one includes the effects induced by non-zero electroweak 



couplings, additional contributions are possible. 

3.2 Corrections to Electroweak Vertices in the MIA 

Now let us consider the effect of supersymmetric contributions to the various electroweak 



vertices in the MIA. As stated in section 2.2, the CKM matrix that appears, for example, 
in the chargino vertex ( |2.33| )-( ^.35| ) is related to the physical CKM matrix by the rela- 
tion ( 2.24| ). At first order in the MIA, the vertex and self energy corrections arising from 
gluino exchange cancel due to SU(2)i x U(l)y gauge symmetry. The first corrections to K 
therefore appear at second order, through diagrams involving two SU(2)l x U(l)y breaking 
insertions on one of the squark lines. The contributions to the vertex therefore tend to be 
suppressed by factors of either m^, or cot (3 and, whilst we take into account these effects 
in our numerical analysis, to a good approximation we may set K = K'^^ . A similar result 



holds for the effective right handed coupling of the W boson (|2.23|) . 

Turning to the Z boson vertex, once again we find that, to first order in the MIA, the 
self energy and vertex corrections cancel due to 'S\]{2)i x U(l)y gauge symmetry The first 
non-zero contribution arises from the diagram shown in Fig. |^ involving two SU(2) l x U(l)y 
breaking insertions. Evaluating the contributions to the effective vertex we find for {C^)^^ 



(^^)23 



92 



asC2 (3) fz {x^) 



2 cos I 



2tt 



ni~ 



[1 + es tan/3) 



{md,LR) 



rribfj, tan (3 
1 + ea tan (3 



Slr 



23 

(3.13) 



where the function fz is given in appendix A^. The expression in square brackets in 
the above expression represents the effect of the flavour diagonal RL insertion. The off- 
diagonal elements of the bare mass matrix can also induce terms proportional to and 
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H- 



H- 



Figure 5: The dominant GFM contributions to the left-handed charged Higgs vertex, arising from 
gluino and higgsino exchange. 



Sjij^, that can viably compete with the corresponding contributions that arise at third order 
in the MIA. Although the vertex ( |3.13| ) is enhanced by tan/?, we shall see later that the 
contributions to a given process due to this vertex scale as 'm'z/^susY typically 
rather small. 

Now let us turn to the Higgs sector, where the effects induced by supersymmetric 



contributions to the charged and neutral Higgs couplings are known to be large |11, 12, 24 1. 
These corrections can, in turn, affect FCNC processes especially in regions of parameter 
space where FCNC mediated solely by sparticle exchange are suppressed by large sparticle 
masses. 

The charged Higgs vertex receives corrections p5| , pi] , [T2| from both gluino and higgsino 
exchange. To second order in the MIA, the effective charged Higgs coupling is given by 



cr 



92 



C 



H+ 
R 



92 



.. cot (3 \k,j (l - tan/3 + e'y(^b°V'53j tan/?) + Af. 



{rrid) 



n 



^mvK (1 + eitan/3) 



(3.14) 
(3.15) 



where i,j = 1,2,3 and vj^^^ = {Y^^')'^:^. The factors e'^ and e'y are given by 



.(0) 



-^C2(3)-^F2 



1 (^d,LR)33 ^ 



2 [yuL^Vd, 



(3.16) 



The arguments of the loop functions H2 can be obtained by the appropriate generalisations 
of ( |3.5D . Finally, the 3x3 matrices A^'^ denote the additional off-diagonal contributions 
that arise in both MFV and GFM models due to the off-diagonal elements of the bare 



mass matrix and the GFM parameters. A^y may be decomposed as follows 

A L,R _ aL,R L,R 
ij ij ^ hj ■ 



(3.17) 



The MFV contributions A^- to the vertex have been highlighted in 14 ] . In the formalism 
developed in section 0, they arise due to the presence of the bare mass matrix in the 
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neutralino vertex, and have the following form 



e^yey(y/°Vttan/3)^ 
1 + ej tan (3 



=0, 



= (3,1),(3,2), (3.18) 
otherwise. (3.19) 



It should be noted that the additional terms found in |14| for = (1,3), (2,3) do not 
appear as we do not assume that the trilinear soft terms are proportional to the bare 
Yukawa coupling. If one adopts the parameterisation described in appendix where such 
a relation is assumed, it can be shown that one obtains an additional contribution to A^- 
in agreement with [l^ . 

The GFM contributions to A^- arise from the two diagrams shown in Fig. ^. Evaluating 
the contributions yields 



Kii tan P 



^'^^ 1 + 6, tan/? 



+ 



m^efiL6^(n^"0%~^,(l + 63tan/3) 
mfo (1 + ej tan fj) 



Kii tan P 



(i,j) = (3,1),(3,2), 

(3.20) 

(i,j) = (1,3),(2,3), 

(3.21) 



where e'j^j^ and e'j^j^ are 



167r2 r^(0) 



-^2(yii„yjJ. (3.22) 



33 



It should be noted that the third term in ( |3.20| ) is proportional the Yukawa coupling of 
the down or strange quark. We include it however as the factors of cos /? present in the 
denominators of the Yukawa couplings (we remind the reader that ~ rrid/mw cos fi) can 
effectively lead 7/^- to vary as tan^ j3. This term can become important if {5'lr)^2 ~ ^RL 
large O (IQ-^). 

Now consider the right-handed coupling of the charged Higgs. In this case the dom- 
inant corrections to the vertex are due to the self-energy correction S^^. The MFV 
contributions to the vertex are reflected by the appearance of a factor of (1 + ejtan/?) in 
the denominator of ( 3.15| ), in agreement with [p!^ ]. 
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In models with GFM it is possible to generate additional terms of the form 
■ eRRX ^^tanp (m^).. 



hj 



(1 + ejtan/3) {md)jj 

TUT, 



3J 



(>RR 



Orl 



= (3,1),(3,2) 



(3.23) 



7^ 



€LLXy tan/3 . 
dL f^d 



{1 + Cj tan (3) 



LL 



{md) 



LR 



(i,j) = (1,3),(2,3). 



(3.24) 



A particularly interesting consequence of (|]2D is that one can often avoid the factor of 
the strange quark mass, that appears in the right handed vertex ( p. 14 ) when i = 3 and 
j = 2, via flavour violation in either the RL or RR sectors. 

It is apparent from the above expressions that GFM contributions, to both the left and 
right-handed vertices, can play the role of the off-diagonal elements of the CKM matrix. 
The off-diagonal BLO corrections to the charged Higgs vertex can therefore be rather 
large in the GFM scenario. Substantial enhancement or suppression of charged Higgs 
contributions to FCNC are therefore possible, even in the limit where the squarks decouple 
from the theory. In addition, tan/3 enhanced corrections affect the underlying structure 
of the charged Higgs vertex, in both GFM and MFV, via the factors of (1 + ejtan/3) that 
occur in the denominator in (|3.15|) and the corrections e'g and e'y that appear in ( 3.14| ). 



The corrections to the charged Goldstone boson vertex [11| prove to be rather small, 
as the vertex is protected by SU(2) symmetry and the self energy and vertex contributions 
approximately cancel. These cancellations are required as, in a general gauge, the 
corrected Goldstone boson vertex must act to cancel the ^ dependence of the contributions 
originating from W boson exchange. The corrected vertex must therefore, in a similar 
manner to the corrected W boson vertex, be proportional to S\](2)l x U(l)y breaking 
effects, even for GFM. 

Finally, let us consider the corrected neutral Higgs vertices ( 2.31| )-( ^.32| ). The dom- 



inant contributions originate from the self energy corrections 5]^^. To first order in the 
MIA the contributions to the flavour diagonal elements of the effective vertex become 



C 



AO 



«52 



tan 13- 



{md)i 



2mw ' (1 + ej tan /?) ' 
whilst the contributions to the effective and vertices are 



92 



(md) 



+ 



2 cos Pmw (1 + tan /?) 

92 {md)ii 



(cos Q + ej sin a) , 
(sin a — ej cos a) . 



(3.25) 

(3.26) 
(3.27) 



2 cos (3mw (1 + tan (3) 

At third order in the MIA, further corrections proportional to combinations of and 
4r ^'^^ generated in a similar manner to ( ^.6| ) that can lead to large corrections to the 



Yukawa couplings of the first two generations. Full expressions can be found in [26| 
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The off-diagonal elements of the coupling are generated by MFV and GFM contri- 
butions and, in a similar manner to the charged Higgs vertex, it is useful to perform the 
decomposition 

Kh)°™. (3.28) 

where = A^,G^, H^, hP . The dominant MFV contributions to the off-diagonal elements 
of the coupling arise from higgsino exchange and are given by 1^, |T4| , 



/ ij 2mvi/ (1 + ej tan/?) (1 + e,- tan/3) 



The GFM contributions arise primarily from gluino exchange and yield the additional 
contribution 



ij 2my^ 



{l + ei tan/?) (1 + ejtan/?) 



"LL 



+ (l + 6,tan/3) V^^)., + (l + .,.tan/?)2 



(3.30) 



The off-diagonal couplings of the scalar Higgs bosons and may be obtained via the 
simple substitutions 

Cf = i sin (a - /3) cf , cf = i cos (a - /?) cf , (3.31) 

whilst the right handed couplings can be obtained by taking the Hermitian conjugate. Due 
to an accidental cancellation between the self-energy and vertex corrections, the terms 
proportional to and vanish at LO. However, once BLO corrections are taken into 
account, it is possible for these insertions to reappear through their effects on the bare 
mass matrix rri'^^ |17|. 



Once again, it should be noted that, due to SU(2) invariance, the Goldstone boson 
vertex does not receive large corrections even once GFM contributions are taken into ac- 
count. As a result any contributions to the corrected vertex are attributable solely to 
SU(2)i X U(l)y breaking effects and are rather small. 

3.3 Additional Electroweak Effects 

As discussed at the beginning of this section the results presented so far have been derived 
in the limit where the electroweak gauge couplings gi and g2 are set equal to zero. The 
aim of this subsection is to briefly discuss the dominant contributions that arise once one 
proceeds beyond that approximation, and to provide some simple substitutions such that 
these effects can be taken into account. 

First, let us consider the effect of such corrections on the bare mass matrix. One of 
the most important corrections in this case is due to the gaugino-higgsino mixing diagram 
shown in Fig. ^ that arises if the insertion (5]^^ is non-zero |21]. The corrections induced 
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-X- 



/ 



— ' — J r\y>y\r<X. i — ^ 



w- 



Figure 6: Contributions to m^°'' arising from the insertion (5|;^ 



by this diagram may be taken into account by making the following substitution in (3.8) 



(3.32) 



a=l 



where y~ = n^^'d llI'^'^ - ^yLL given by 



47r sin2 Ow V^mw sin (3 \ '^^ ' 



(3.33) 



a denotes the electromagnetic coupling constant. We have made use of the relation (|2l2^ ) 
to express the contribution in terms of flavour violation in the down squark sector. In the 
phenomenologically interesting region /i > and At < 0, interferes destructively with 

the gluino contribution sll and acts to reduce the correction to m^^^ that is proportional 
to flavour violation in the LL sector. 

Turning to the charged Higgs vertex, as discussed in [14|, large contributions to the 
left-handed vertex arise from diagrams featuring gaugino and higgsino exchange. They 
may be included by making the following substitution in ( |3.14 ) 



4 + 4- 



(3.34) 



e' has the following form 



a 



m^o — f 2 



Att sin"* 9w ^ 

^' a,a Xa 



-Mr 



-K*2Kitan%i/2(y2^,0 



+ 



(3.35) 



where the quantity Maa is given by 



Our results for agree with those originally given in |14]. To include the additional effects 
induced by GFM, one has to consider the diagrams shown in Fig. ^. Their effects may be 
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Figure 7: Additional GFM contributions to the left-handed charged Higgs vertex, due to the 



electroweak effects considered in subsection B.S 



included by making the following correction to ( |3.20| )-( |3.21| ) 



(3.36) 



where eZ'n given by 



a 



txLL 



47r sin^ dw 



E 



m - 

a Xa 



N. 



a2 



V* N* 

V2 



(3.37) 



Both the MFV ( p. 35 ) and GFM ( |3.37| ) corrections typically interfere destructively with the 
dominant gluino contributions and can lead to an appreciable reduction of BLO effects. 

Finally, let us consider the neutral Higgs vertex. As discussed in the previous sub- 
section, the dominant contributions to the effective vertex arise from the self-energy cor- 
rections SJ^^. The effects induced by non-zero electroweak couplings may therefore be 
included, in a similar manner to the bare mass matrix, by making the substitution ( 3.32| ) 



in djOj) 



3.4 Other Methods 

The method outlined in section || takes into account both tan j3 enhanced effects and those 
induced by non-minimal sources of flavour violation. Other methods have been proposed 
in the literature that can be modified to include the effects of GFM and it shall be useful 
to briefly consider how two specific examples compare with the method employed in this 
paper. 

The first method, presented by Buras et al. works in the bare SCKM basis. In this 
basis the Yukawa matrices vjf''^ and 1^^'^^ that appear in the superpotential are diagonal. 
Calculating the self-energies in this basis gives the physical quark masses 



(3.38) 



where Y^^^ is the diagonalised Yukawa matrix, and dih^^^ denotes the contributions of the 



self energy corrections (2.5) calculated in the bare SCKM basis. Dl and denote the 
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unitary transformations performed on the squark fields that transform between the bare 



and physical SCKM bases. The bare mass matrix m!"^'^ defined in 
quantities by 



m 



(0) 



A- 



is related to these 



(3.39) 



It is straightforward to relate the matrices to the unitary matrices Vd^ ^ that ap- 

peared in section |2|. If, in analogy with the transformations ( |2.8D -(^^), one defines a 
transformation from the interaction basis to the bare SCKM basis such that 



v(o) 



(0) (v^o)o\ ym 



' Y 



Dr and Dl are then given by 



DL = Vd,V,^ 



(o)t 



One may also define the bare CKM matrix in the SCKM basis 

i^(0) = v;Wyjo)t ^ uIkdl. 



(3.40) 



(3.41) 



(3.42) 



The role played by the off-diagonal elements of rn'^'' in section |2| is taken by the unitary 
matrices D^, Dr and the bare CKM matrix K^^\ For instance, the bare CKM matrix 
elements and , in the bare SCKM basis, are related to the corresponding matrix 
elements K, defined in the physical SCKM basis via the relation 



(0) 
ts 



(1 + £3 tan /3) 
"(1 + e, tan/3) 



ts 



(1 + e3 tan/3) rUg _^ ,d _ ^ll tan /? . 
'"^^ (1 + e, tan (3) rub ''~^rl^^^ (1 + tan /3) ^^^^^^ 



K 



(0) _ (1 + estan/g) 
(1 + es tan/3)' 



cfe 



(1 + £3 tan/?) rrig 
^RL — — X 



[\ + estan/3) 



-^-^^^^ (l + e.tan/3)''^^'^^^ 



Ktb, 
(3.43) 

(3.44) 



Where we have used the shorthand 5^^. 



Kl)23 and 



('^Lij)23' Strictly speaking, 
the uncorrected CKM matrix K should appear in the above relations, however, as discussed 



in section |3.2| , the vertex and self energy corrections are negligible and one may, to a good 
approximation, set K = K'^^. An interesting consequence of this formula is that the matrix 



element K^^^ obtained by diagonalising the bare Yukawa couplings Y^"^^" can be zero in 



(0)o 



the presence of general flavour mixing |21, 16|. We will discuss the consequences of this in 



section 7.4 



As the two methods are practically equivalent, choosing between them essentially be- 
comes a choice as to which is more suitable for the problem at hand. In MFV scenarios the 



method presented in [14| is generally more convenient as it is only necessary to calculate the 
diagonal parts to the vertex and self-energy corrections induced by gluino exchange. For 
example, when using the method described in section || the correct form of ( |3.29| ) is only 
obtained when one considers the off-diagonal gluino contributions as well as the higgsino 
exchange diagram. 
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In the GFM scenario the situation is rather different. As the off-diagonal contributions 
to the electroweak vertices and non-renormahzable operators induced by gluino exchange 
are evaluated anyway, the method described in this paper can become more preferable. 
In particular, the flavour diagonal contributions, induced by the exchange of the super- 
symmetric particles, to the various non-renormalizable operators applicable to the process 
under investigation no longer have to be calculated, as the role played by the matrices Dl 
and Dr is replaced by the off-diagonal elements of m^^^ . 

The second method, presented by Dedes and Pilaftsis |jl^ , concerns itself mainly with 
CP violation, however it is essentially applicable to both MFV and GFM CP conserving 
scenarios as well. After translating to the physical SCKM basis their expression for the 
bare mass matrix reads 

ma = mf^R, (3.45) 

where is a 3 x 3 matrix. It is then possible to express the various Higgs interactions via 
an effective Lagrangian expressed in terms of the physical quark masses and the inverse 
of R. In the context of MFV and GFM scenarios with mixing only in the LL sector this 
parameterisation is sufficient. However taking into account all sources of flavour violation 
yields the more general form 

rud = mfRL + RRTuf + Rg- (3.46) 
In the MIA, the 3x3 matrices Rl, Rr and Rq may be decomposed in the following way 
R^ = 1 + est&n [3 + ^eYK^Y!^K tan 13 + eLLtan fi x^ji^, 

u 

Rr = +eRR tan /3 x^jj^R, (3.47) 
Rg = -eRL rUgX^^jjiL- 

Obtaining a solution for m^^'^ is therefore rather more complicated than simply finding the 
inverse of R. Considering each element of m^^^ in turn, however, it is possible to replicate 



the results for m^^^ presented in subsection ^A. 



4. B ^ Beyond the LO 

Of all the FCNC processes involving transitions between the b and s quarks, B Xgj is 
currently the best understood both experimentally and theoretically. The data being taken 
by B-factories such as BaBar and BELLE, is leading to an increasing degree of precision 



for the measurement of the branching ratio of the decay. The current world average is |27| 



BR {B ^ X,7)^^p = (3.39;°;30) x IQ-^ (4.1) 



This value takes into account the most recent BELLE [^] and BaBar results [2S]. 
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The SM prediction for the branching ratio is based on a NLO calculation that was 



completed in Refs. |3^, 31, resulting in the prediction 



BR {B ^s7)gM = (3.70 ± 0.30) x lO""^. (4.2) 

It has been pointed out recently [^] that, if one applies a realistic cut-off for the 
photon energy (rather than E-y > 1/20 m^), a dependence on two additional energy scales 
appears when calculating the branching ratio for the decay. The first scale (//j = ^/mf,A) is 
associated with the energy of the final hadronic state Xg , whilst the second is dependent on 
the energy range under investigation (^uq = A = mf, — 2E^). The perturbative uncertainties 
associated with these scales are rather large and can lead to a significant increase in the 
error associated with the branching ratio. However, in exchange, the final result can 
be compared directly with those determined directly from experiment rather than model 
dependent extrapolations to Ej > 1/20 mi,. 

We should also briefly mention that steps are now being taken towards a NNLO cal- 
culation 1 34, 35 1, that should increase the accuracy of the SM prediction to roughly 5%. 

The effective Hamiltonian relevant to AF = 1 processes such as the decay B Xg'j 

is 

Heff = -^K!^*K!^ ^ [a (m) if,) + C[ (/x) 0[ (/x)] . (4.3) 
The operators most relevant to the decay are 

= -^rn, {-SLa^^^hR) F^,, = {sLa''''T%n) G%. (4.4) 



(The six remaining operators can be found, for example, in [£0|). The primed operators 
can be obtained via the simple substitution L <-> i?. The contributions to the primed 
operators are negligible in the SM. However, in more general models, such as the MSSM 
with general flavour mixing, their effects can no longer be ignored. As the primed and 
unprimed operators do not interfere with one another, any new physics contributions to Cy 
and Cg enter quadratically and therefore act to increase the value of the branching ratio. 
New physics contributions to C7 and Cg, on the other hand, interfere directly with the SM 
contribution and can lead to far more varied effects. 

The good agreement (within \a) of the SM prediction and the current experimental 
results allows one to place increasingly stringent bounds on the effects and mass scale of 
new physics contributions. In doing so it is important to include the effects of new physics 
at a similar precision to the SM result. NLO matching conditions have been completed 
for several extensions of the SM, for example, the NLO matching conditions relevant to 



the 2HDM were presented in |37, 36 1 whilst a more general analysis was presented in 
Turning to the MSSM, however, things become rather more complicated. A complete 
NLO calculation would involve the evaluation two loop diagrams involving both gluons 



''This result includes the NNLO effect induced by using the running charm quark mass rather than the 
pole mass when calculating the charm quark contributions to the decay . A more formal NNLO analysis 
of these effects has been performed in P5| 
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and gluinos. For MFV this task is already underway and, for example, the NLO matching 



conditions for the charged Higgs contribution have been discussed in [39|. Theoretical 
calculations have, thus far, concentrated on particular cases. The calculation presented 
in 1^^, for example, considers a realistic but specific region of MSSM parameter space 
where the charginos and lightest stop are relatively light compared to the rest of the 
sparticle spectrum and tan [3 is rather small. These results were extended to the large 



tan/3 regime in [11, 12 1. The same papers also considered the dominant effects that occur 
BLO for generic SUSY scenarios taking into account effects enhanced by large logs and 
tan/3. These results were subsequently extended to include CP violation [^], additional 
CKM effects Q and SU(2)l x U(l)y breaking and electroweak effects in [||]. 



In the GFM scenario the LO matching conditions have been known for some time |42 



H, H] and, in a similar manner to the MFV calculation, the NLO matching conditions 



have been derived in the limit where the gluino decouples and tan/5 is small [38|. An 



extension of the analysis given in [11| to the GFM scenario was presented in 16 1 where 
it was found that BLO corrections can play a large role and can lead to a significant 
relaxation of the limits placed on GFM parameters compared to a LO analysis. The aim 
of this section is to present calculations in the MIA for both the electroweak and SUSY 
contributions to B — > X<j7, allowing one to easily determine the dominant effects that occur 
once GFM is taken into account compared to MFV calculations, as well as presenting the 



calculations detailed in |16| in a more transparent way. In doing so we therefore adopt the 



approximations discussed at the beginning of section y. 

4.1 BLO Corrections to Electroweak Contributions in the MIA 

As emphasised in section ^, including tan jS enhanced corrections to the charged and neu- 
tral Higgs vertices can lead to large effects. Let us first consider how the LO charged 
Higgs contributions to the decay B — > Xg^ are altered once these effects are taken into 
consideration. Using the corrected vertices ( |3.14D ~( ^.15 ) the dominant MFV contribution 



in the large tan f3 regime is given by |11, ^, |I 



' / l + e3tan/3l " ^ l + e,tan/3 i ^'^ V"^^ 



The loop function Fj^^ is given in appendix A. 2. Note that ( [4.5|) includes the LO contrl 



bution in addition to the BLO corrections. Turning to the degenerate mass limit we see 
that depends on the sign of //. In the phenomenologically favoured region ^ > 0, for 
example, the BLO corrections induced by gluino exchange typically reduce the branching 
ratio compared to a simple LO calculation. The higher order contribution proportional to 
Y^l^'^'^ , first pointed out in on the other hand is dependent on the sign of the trilinear 
soft terms and can therefore interfere destructively or constructively with the (dominant) 
gluino correction depending on the model at hand. It should be noted that ( [4 .51 ) only 
serves as a rough approximation of the BLO effects and that the additional effects arising, 
for example, from gaugino mediated exchange and SU(2)i x U(l)y breaking can lead to 
deviations from this idealised result in some regions of parameter space |14]. The dominant 
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effects that arise from these corrections may be included by performing the substitutions 



presented in subsection 3.3 



The GFM contributions to the charged Higgs vertex, discussed in section 3.2, give rise 
to additional BLO corrections to Cr^s given by^ 



5" C7,s 



GFM 



tb 



ts 



^LL^d, ^'^^ ^Y^LLX^ {Y^°^ tan /3) 



(l + estan^) (1 + tan /3) (1 + es tan/3) 



+ 



nib 



1 + tan j3 



LR 



y^yW tan /3 4L kS 



K<1 

(>LL 



m 



(4.6) 



The factor of K^f^/Kf^ that appears in front of ( [4.6[ ) reflects the fact that the flavour 
change is governed by the GFM parameters S'li^ and 5^^, rather than the CKM matrix. 
The additional GFM contributions interfere directly with the MFV corrections to the LO 
result and, depending on the sign of 6f^^ or Sf^^, can easily lead to large reductions or 
enhancements of the MFV result. In addition to these contributions to Cy^g it is also 
possible, for non-zero (5^^ and to induce corrections to the primed Wilson coefficients 



GFM 



mt 



2-^7,8 



mt 



m. 



, [l-(4 + 6V(yr)2)tan/3] 



K* 



(l + es tan/3) 
(l + es tan/3) 



(2) 



7,8 



m-p, 



mt 



m 



H+ 



(4.7) 



As the LO contributions to the primed coefficients are suppressed by factors of mg/mi, the 
dominant behaviour, once BLO corrections are taken into account, is determined solely by 
GFM effects. Note that these GFM effects persist even if the squarks decouple from the 
theory. 

It has been pointed out in Refs. |3|, |l^ that the corrected neutral Higgs vertex can 
also induce contributions to Cj^s through the diagram where a neutral Higgs boson and 
a bottom quark undergoing a chirality flip are exchanged. In the limit of MFV using the 
corrected vertex ( |3.29| ) one obtains ||3|, |l^] 



S""c,,s) 



= - 14 tan' /5 

36 mj^ (1 + estan/3) (1 + e3tan/3)'' 



(4. 



The tan^ /3 dependence of the Wilson coefficient is characteristic of the corrected Higgs 
vertex ( p. 291 ) and can compensate for the suppression factor m^/m^. 



''From now on we shall adopt the conventional shorthand Sxy ~ {^xy) ^ 
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Turning to the effects of GFM contributions, using ( 0.30 ) it is possible to induce 
additional corrections to Cr^s 



/ 36 m^KLK^h 



GFM 1 rn 



dL 

(1 + estan/3) (1 + e3tan/3)^ 



rrif, (1 + tan/?) (1 + £3 tan /?) 



(4.9) 



In a similar manner to (^^), the GFM contributions arising from neutral Higgs exchange 

vary as tan^ f3. The primed coefficients also receive contributions if 5^^, ^'rb. 

zero 



rn~g ^Riesx^^ . 

n oOpr 

(1 + 63 tan /?)2 

4.2 BLO Corrections to SUSY Contributions in the MIA 



(4.10) 



The super symmetric contributions to the decay B X^'j can proceed through a number 
of channels. In MFV, the only SUSY contributions arise from diagrams involving chargino 
exchange. Once GFM effects are taken into account, additional diagrams arising from 
FCNC mediated by gluinos and neutralinos can occur and give rise to contributions to both 
the unprimed and primed Wilson coefficients. As the gluino contributions are enhanced by 
factors of q<j (compared to the MFV contributions), these effects are rather large and can 
play an important role for even small deviations from MFV. 

All four insertions give rise to contributions to either Cy^g or their primed counter- 
parts and it will be useful, for our purposes, to decompose the overall gluino mediated 
contribution to the decay as follows 

/ N MFV nLL nLR nRL nRR 

53^7,8 = {S9Cr,s) + (S'Cr^s) + [S'Cr^s) + {s'Cr,s) + [S'C-r^s) (4.11) 

The primed coefficients and other SUSY contributions may be defined in a similar manner. 
The dominant BLO corrections to the gluino contributions, shown in Fig. ^, arise from the 
flavour violation mediated by the off-diagonal elements of the bare mass matrix and are 
proportional to m^^^ tan (3. 

The MFV terms present in the bare mass matrix ( |3.8| ) can lead to a correction to the 
gluino contribution of the following form 

,~9r \ _ 8 a, sin^ Owml, eyY^ tan^ /? /x [7,8] / n 

^^'^V - 3a m| (l + e,tan/3)(l + e3tan/?)m5'6 ^'^'^^^ 

[7 01 [7 01 

The loop functions 1\ ' (x) and Jj ' (x), that appear throughout this section, can be 
found in appendix A. 2, denotes the ratio 

2 

m~ 

x~^ = ^ (4.13) 
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Figure 8: BLO corrections to C7 (Cs) and C7 (Cg) arising from gluino exchange, the photon 
(gluon) line is attached in every possible manner. 



where m- is a common mass of the quadratic soft terms ( p.l2D (that is. 



m| 



m 



d,LL , 



) ) . Due to the 1 / m~ suppression of the amplitude it would be expected that these 
effects are typically rather small when compared to the additional BLO effects arising, for 
example, from the modified charged Higgs vertex. However, the terms feature a dependence 
on tan^ (3 in the numerator and should be included if one wants to consider the effects of 
all tan/3 enhanced corrections. We should note that this correction is entirely consistent 
with the definition of MFV presented in Q and is a result of the transition between the 
bare and physical super-CKM bases. 

The GFM contributions to the Wilson coefficient Cj^s arising from gluino exchange are 
due mainly to the LL and the LR insertions. The contributions due to the RL and the RR 
insertions are suppressed by factors of the strange quark mass and may be safely ignored. 

Contributions arising from the insertion are generated at first and second order in 
the MIA. At first order, the chirality flip is generated via the bottom quark that appears in 
the operators ( |4.4| ). At second order, the contribution arises from the diagram involving a 
diagonal LR insertion, an LL insertion and a chirality flip on the gluino line. This correction 
can play an important role for even moderate tan /?, and for large tan /? dominates the 
overall behaviour of the contribution to C7. If we ignore the effects generated by the 
diagonal elements of the trilinear soft terms, we have 



LL 



9 / \ 2 

8 Os sin 6w I my/ \ 



a 



m~g J 



CLL tan^ P 



(1 + eg tan (3) (1 + €3 tan /3) rrig 

•mill eiL tan^ /3 

mi (1 +e,tan/?) (1 +e3tan^) 



-'6 



i-d) 



+ 



tan 13 



y^d) 



2mg (1 + es tan (3) 



(4.14) 



The first and second terms in square brackets that appear in ( |4.14 ) arise at the respective 
order in the MIA. The chirally enhanced BLO term (that is proportional to 1^'^'^ {-^d)) 
occurs at first order in the MIA, due to the off-diagonal elements of m^^^ . This term tends 
to reduce the overall effect of the contribution that arises at second order in the MIA (the 
term proportional to Jg^'^^ (-^d)) ^ > (this is one of the contributions to the focusing 
effect discussed in |15, |l6[). For ^ < on the other hand, the two contributions interfere 
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constructively and increase the contribution to relative to a LO calculation. The LO 
contribution proportional to I5 (x) also undergoes a similar reduction once BLO effects 
are taken into account. However, in this case, the BLO correction is reduced by a factor 
of ml/m~. 

For non-zero Sf^p,, the dominant contribution at LO arises from the diagram involving 
an LR insertion and a chiral flip on the gluino line. This contribution is therefore en- 
hanced by a factor of m-g/mi,- Higher order contributions in the MIA do not feature this 
enhancement and are typically rather small. To second order in the MIA we have 



LR 



+ 



8 sin 6w Os 
~ 3K*^Ktb ^ 



tan P 



nig 



1 



2m| (1 + tan ^) (1 + £3 tan /?) 



rub (1 + eg tan P) 



[7,8] 



LR 



(4.15) 



Once again the first and second terms in the square bracket arise at the respective order in 
the MIA. From ( p^ it can be seen that BLO effects can reduce, or enhance, the dominant 
contribution due to the insertion that arises at first order in the MIA, depending on 
the sign of e^. In the phenomenologically favoured scenario ;U > 0, in particular, is 
positive and BLO effects act to reduce the LO contribution to Cj^g. The term that occurs 
at second order in the MIA tends to be subdominant, compared with the chirally enhanced 
term that appears at first order, but acts to reduce the contribution to Cy^g further. 

Turning to the primed sector, the corrections due to MFV, LL and LR contributions 
are suppressed by factors of rris and are rather small. We are therefore left with the 
contributions arising from RL and RR insertions. 

The contribution due to the insertion Sf^^ to second order in the MIA is given by 



RL 



8 sin 0w Os 



mw\' 
9 / 



m 



tan/3) [7^; 



nih (1 + estan/?) 



i-d) 



fJ-rrib tan/3(l +eyyf^tan/3) [7^8] , 



2m| 



(l + egtan/?)^ 



(4.16) 



Comparing the above expression with (4.15) we can see that the form of the two are rather 
similar, the only differences being the replacement of with £3 in the denominator of 
(4.15) and multiplication by an overall factor of 1 + eyY-l^ tan p. We therefore see that BLO 
corrections, once again, act to reduce the contribution due to (5^j^ with respect to a purely 
LO calculation if > 0. 

Finally the contribution to Cj g arising from non-zero has the form 



\RR 8sin2 6'H^a, 



tan P 



2mg (1 + 63 tan/3) 



J. 



[7,8] 



i-d) 



r[7,8] 



i-d) 



eRRt&n^ P fi [7,8] 



mlfi 



(l + e3 tan Pfrug 
eRR tan^ /? 



i-d) 



mi (1 + 63 tan/?) 



(x- 



(4.17) 
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Figure 9: BLO corrections to C7 (Cs) and C7 (Cg) arising from cliargino exchange, the plioton 
(gluon) line is attached in every possible manner. 



In a similar manner to ( [4.14| ), the chirally enhanced BLO term arising at first order in the 
MIA, due to the ofF-diagonal elements of the bare mass matrix, can once again affect the 
dominant, chirally enhanced, LO contribution that arises at second order in the MIA (the 
first and second order terms that appear in the square brackets respectively). 

We now turn to the chargino contributions to Cj^s and C^g, and use an analogous 
decomposition to ( 4.11| ). The dominant BLO corrections to chargino exchange arise from 
the diagrams shown in Fig. ^. 

The effect of BLO corrections in the limit of MFV are well known |11, 14|. In GFM 
models, however, it is possible that additional sources of fiavour violation in both the up 
and down squark sectors can significantly alter the MFV result. 

Contributions from flavour violation in the up squark sector enter at LO and are 
therefore rather large. As we are chiefly concerned with flavour violation in the down 
squark sector, the only relevant source of flavour violation in the up squark sector is the 
insertion (5^^, which is related, by SU (2) symmetry, to ( |2.2[1| ). Flavour violation in 
the down squark sector, on the other hand, only enters via BLO effects induced by the 
off-diagonal elements of m^^'^ . 

The contributions that arise due to LL insertions have the form 



6^ O 



7,i 



m 



Xa 



LL 



K* 
K* 



cs_ \ " 

'* / -I 
a=l 



m 



w 



m 



Va 



Xa 



mw \/2 cos /3 (1 + es tan [3) 



-'2 



LL 



+ 



Kcs mw 



U*2VaieLLX^tanP 



Ktb m^- ^2 cos Pil + €s tan/3) (1 + £3 tan f3) 



H. 



[7,8] 



"LL 



(4.18) 
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Once again the loop functions H'f'^'^ and /j"^'^' can be found in appendix [A.2| , the ratio y~ 
is given by 



2 

m~ 



m 



(4.19) 



Xa 



The dominant contribution at LO for large tan (5 arises from the second term in square 
brackets in ( [4.18| ) that is enhanced by factors of m^- /m\\r and 1/ cos (3. This term, for > 



0, has the same sign as the gluino contribution ( [1.14 ) and acts to enhance the contributions 
due to flavour violation in the LL sector. The BLO corrections to the Wilson coefficient 
are reflected by a factor of (1 + estan/?) that appears in the denominator of the second 
term in square brackets in (|4.18 ), and the term that appears in the second line of ( 4.1^ ), 
proportional to Both of these BLO corrections for > act to decrease the LO 

contribution. For /i < 0, on the other hand, both corrections act to increase the Wilson 
coefficient relative to the LO result. 

The LR insertions also contribute to the Wilson coefficient 



LR 



m~g 



eRLX^ 



KtbKl ^ ^/2m^~ cos (3 m,, ( 1 + e, tan /3) 



X H. 



[7,i 



(4.20) 



In this case, GFM contributions only enter via BLO effects induced by the off-diagonal 
elements of the bare mass matrix. For /i > and large tan/3, the contribution ( 4.201) has 
the opposite sign to the gluino contribution ( [4. 151 ) and large cancellations are possible, 
which contribute significantly to the focusing effect pointed out in [^]. 

Before proceeding with the results relevant to the primed sector we should note that 
once again, the contributions to Cy^g arising from RL and RR insertions are suppressed by 
factors of mg. 

In a similar manner to the corrections that arise from gluino exchange, the only dom- 
inant contributions to the primed coefficients arise from RL and RR insertions. The con- 
tribution arising from RL insertions is given by 



6^ C7' 



RL 



E 



K* 



1 ^ts^Xa 



- — 2~'^a2'^a2tan/?i2 

2m „ 



Xa 



mw 
V2cosP 



V:^Ua2 ( H. 



r[7,8] 



mb (1 + estan/?) 



(4.21) 



and the contribution arising from RR insertions is 



RR 



E 



K* 



a=l 



Kim 



Xa 



;r^^a2f^a2tan^/^ 
2m „ 

Xa 



mw 



V2COS/3 



eRRXit&nP 



(1 -hes tan/3) 



(4.22) 
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LO contributions to both coefficients arising from either MFV or non-zero SJ^j^ are sup- 
pressed by factors of ms/nib and are therefore rather smah. At large tan/3, therefore, the 
BLO effects dominate the behaviour of the chargino contributions to the primed opera- 
tors. We should also note that in a similar manner to the case of the LR insertion both 
corrections (for fi > 0) have the opposite sign to the gluino contributions ( 4.16D -(4.17). 

In the GFM scenario, neutralino contributions also play a role. They can become 
especially important when, for example, the gluino and chargino contributions partially 
cancel [U]. However, as the neutralino couplings (|C.6D -( |C77|) are rather complicated, we 
shall refrain from presenting complete analytic expressions for the coefficients in the MIA. 
The overall effect of including BLO corrections is to modify the Wilson coefficient in a 
similar manner to the gluino contributions ( 4.15| )~ (|4.17| ). As an example, the contribution 
arising from LR insertions to Cj^s due to bino exchange becomes 

The loop function I^^'*' can be found in appendix [A. 2 whilst its argument is given by 
= m-ZMf. The contributions induced by neutral gaugino-higgsino mixing, on the 
other, hand are more complicated due to the appearance of the bare quark mass matrix in 
the couplings (|a6|)-([a7D. 

Finally, let us consider the effect of evolving these coefficients from the SUSY matching 
scale fJ-suSY to the electroweak scale iJ,w- Considering only the mixing between and Og 
we have the LO relation |47| 

(fiw) = v^^Cr {^J'SUSY) + - itJ'SUSv) ■ (4.24) 

The factors of rj in the above expression reflect the resummation of leading logarithms and 
are given by 

77 = — , (4.25) 

where (fisusv) (fJ-w) should be evaluated with the QCD /? function relevant for 

six active flavours. If we retain only the flrst logarithm that appears when expanding the 
factors T] we have 



4 

6C7 (nw) = SC-j ifisusv) - i^as ifiw) 



^Cr ifJ-susv) - g'^^s {^J■susY) 



log^%F 



(4.26) 



From the above expression, it is apparent that the evolution of the coefficients from fisusY 
to iiw acts to reduce the overall SUSY contribution. In addition mixing with the coefficient 
Cg can also play a role. If 6Cs has the opposite sign compared to 6Cj, for example, further 
reductions are possible. 

Finally, let us provide a recipe for implementing BLO corrections into existing LO 
gluino matching conditions calculated in the MIA When performing LO calcula- 

tions, one ignores the corrections to the bare mass matrix, discussed in section 2.1, and 
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the F-terms are therefore assumed to be flavour diagonal. Once one proceeds beyond 
the LO however, this assumption no longer holds and additional off-diagonal elements in 
the squark mass matrix appear, that are attributable to the factors of rri'^^ that appear 



in ( p. 14 ). In the LR sector, in particular, the off-diagonal elements of the bare mass matrix 



are enhanced by a factor of tan (3. The effect of these BLO corrections can be included in 
existing LO expressions by making the following substitutions 

/itan/3(my^ty 

^in) . . - Ur) . . - : , , (4.27) 



n 



/xtan (3 (m^^^ 

. . - . . - . . '\ ■ (4.28) 



^Irr).. m,LL 



n 



Similar substitutions exist for the insertions and (5^^, however the effects are typically 
proportional to C("i^/?7i|,^^y) and may therefore be safely omitted. In each of the sub- 
stitutions given above, the off-diagonal elements of the bare mass matrix are enhanced 
by a factor of tan (3 and can therefore play a rather large role. One should also note that 
the factor of the down quark mass, that appears in flavour-diagonal LR mixings, should 
also be replaced by the appropriate element of m^^\ Following this recipe, it is relatively 



easy to modify the LO calculation presented in, for example, |45] and obtain results in 



agreement with those presented above. We should note that, provided one calculates m^^^ 
to a similar precision, the substitutions can be used to all orders in the MIA. One may 
also use these substitutions in LO expressions for the chargino and neutralino matching 
conditions, however here one must also take into account the factors of the bare mass ma- 
trix that appear in the chargino and neutralino vertices. Finally, let us emphasise that the 
substitutions (l4.27| )- (|4.2g| ) do not amount to a redefinition of the 5's given in (|2.18|) - (|2.19|) 



but merely represent the form of the BLO corrections to LO expressions. 
4.3 Full Calculation 

With our results derived in the MIA in mind let us now outline the steps required to im- 
plement BLO corrections in the general framework outlined in section ^ where the squark 
mass matrices are diagonalised numerically. After performing the iterative procedure de- 
scribed in subsection p.4| one may obtain the BLO charged Higgs and SM contributions 
by using the matching conditions presented in [37, p6|, to account for the NLO gluon 



contribution, and using the corrected vertices presented in section |2| to evaluate the LO 
matching conditions. As discussed in subsection |3.2| the corrections to the SM contribu- 
tions tend to be rather small and can generally be ignored. The effect of the neutral Higgs 
contribution can also be included by using the matching conditions presented in |14|. 

Turning to the super symmetric contributions, BLO corrections can be incorporated 
by using the super symmetric vertices detailed in section 2^ and appendix ^ in the LO 



matching conditions given in [38, It should be noted that one should use the uni- 



tary matrices T^i ^ that are obtained at the end of the iterative procedure described in 
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subsection 2A when evaluating these contributions. One may also, with care, include the 
additional NLO effects that appear if the gluino decouples by using the matching condi- 
tions presented in |^8|. Once one has evaluated all the supersymmetric corrections they 
may be evolved from the SUSY matching scale to the electroweak matching scale using the 
NLO six flavour anomalous dimension matrix presented in |47|. 

With the supersymmetric and electroweak contributions evaluated at the scale fiw it 
is finally possible to calculate the branching ratio for B — > Xg^ according to |^0|. Let us 
note that this recipe is quite general and may be applied to any other process providing 
the relevant matching conditions and anomalous dimension matrices are available. 

In summary, in this section we have included all the relevant corrections that appear 
beyond the LO in the MIA. In the electroweak sector, we have seen that the additional 
GFM contributions to the charged and neutral Higgs vertices can lead to potentially large 
modifications to the MFV results depending on the sign of the insertion at hand. The 
interplay and cancellations between the various supersymmetric contributions, has also 
been shown to be significant once one proceeds BLO |15| and leads to a focusing effect 
in the phenomenologically interesting region > and At < 0. For the insertions 
(5^^ and (5^^, in particular, large cancellations can arise between the gluino and chargino 
contributions to the decay. For the insertion on the other hand, the cancellations play 
a more minor role, as a LO correction to the chargino correction already exists (arising 
from the insertion S'^j^) and tends to reduce the effect of BLO corrections. Finally we 
have seen that the RG evolution of these corrections can lead to further reductions to the 
supersymmetric contributions to the decay. 



5. Bs — > M"*"/^ Beyond the LO 

As B-factories do not run at the energy required to produce large quantities of Bg mesons, 
the best experimental constraint on the rare decay Bs — > is provided by collider 

experiments. The current 95% confidence limits provided by the CDF and D0 |^9| 
experiments at the Tevatron are^ 

BR {Bs ^ Ai+^")cDF < 2.0 X 10^^ (5.1) 
BR {Bs Ai^^")D0 < 3-8 X 10~\ (5.2) 

CDF and D0 intend to further probe regions of up to O (10-^). At the LHC, on the other 
hand, branching ratios of up to O (lO~^) are easily obtainable after a few years of running 
at ATLAS, CMS and LHCb |2|. 

Theoretically, the decay Bs provides one of the cleanest FCNC AF = 1 decay 



channels. It is described by the effective Hamiltonian |22 



-Heff = -^Ktf*KfY. [c^i^^)o^{^A + c[{M{f^)] (5.3) 



* These results are preliminary, one can find the most recent published results in 
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where the operators Oj are given by 

Os = rrib {SaPRba) ill) , O's = ms (SaPLba) {II) , 

Op = rrib {s^PRb^) {h^l) , 0'p = ms {sMa) {Hi) ■ (5.4) 

As the anomalous dimensions of all three operators are zero, the RG running is trivial and 
the overall branching ratio for the process l = fi is given by 

BR {B. ,^,-) \K!!*K!^\^,/T^^ 

X [(1 - Ami) \^s\^ + 1^^ + 2"i'^iol'] , (5-5) 
where rh^ = m^/mB^ and the dimensionless quantities Fi are given by 

Cs,pmb - Cg pTUs 



Fs,p = ruB, 



nib + ms 



FiQ — Cio — CiQ. 



It should be noted from the above expression that the Wilson coefficient of the operator 
Oio is helicity suppressed by a factor of as the Bg meson has spin zero. The SM 
contributions are only proportional to Oio as the Higgs mediated contributions to Os,p 
can be safely neglected. The SM contributions to Cio have been evaluated to NLO |53| 
resulting in the branching ratio 

BR {Bs Ai+^")sM = (3-46 ± 1.5) x 10"^ (5.6) 

The large uncertainty is mainly attributable to the hadronic matrix element /b^ that can 
be determined from either lattice or QCD sum rule calculations. A representative value 
for fB, is^ 

fBs = (230 ± 30) MeV. (5.7) 

In scenarios beyond the SM, particularly SUSY with large tan /?, the contributions 
to Cs and Cp arising from neutral Higgs penguins can become large and dominate the 
underlying behaviour of the branching ratio. Studies in the MSSM have focussed on both 



MFV ||5^, |6y] and GFM scenarios |6^, 17 1 where the corrections induced by the 
corrected neutral Higgs vertex ( ^.29D -( ^.30| ) lead the branching ratio for the decay to 
vary as tan^ /? (for a review see |2j]). At large tan/3 it is therefore quite possible for 
BR {Bs fi'^fi~) to be enhanced by a few a orders of magnitude compared to the SM 
value, providing an ideal signal for physics beyond the SM. The aim of this section is to 
present analytic expressions for the contributions to Cs and Cp in the MIA that include the 
BLO effects discussed in section |2|. We also discuss briefly the effect BLO contributions have 
on the subdominant contributions that arise from box diagrams mediated by neutralinos 
and charginos. Finally we discuss the application of the recipe given at the end of the 
previous section to Bs A*"*"/^"- 

^The current unquenched lattice calculations for fs^ vary from 215 MeV up to 260 MeV (for 
details of the errors associated with these values we refer the reader to the original papers). As the branching 
ratio is proportional to the square of /b, we decide to take a rather conservative estimate for the magnitude 
of as recommended in . 
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Figure 10: The neutral Higgs penguin contributions to Cs,p (on the left) and C'g p (on the right). 



5.1 BLO Corrections to Electroweak Contributions in the MIA 

As above for B Xg'j, here we will present MIA calculations for the contributions that 



arise due to the effective vertices presented in section . 

Corrections to the effective Z vertex ^.13 lead to contributions to Cio proportional to 



GFM 



as fJ-rrib 



tan P 



2>KtbK*^ a mi {1 + e, tan /?) (1 + €3 tan /3) 



fz {x^) Sin (5.8) 



whilst there is a similar contribution proportional to Sj^^ to C[q. Terms proportional to 
may be generated at third order in the MIA, which undergo BLO corrections from terms 
that appear when using the bare mass matrix (|3.8|). As stated in the previous section 



however, Cio and C[q are both helicity suppressed by factors of and their contribution 
to the overall branching ratio is therefore typically limited to the low tan /? regime. 

The charged Higgs contributions to Cio and C(q arise from Z penguins and box dia- 
grams. The contribution to Cio is suppressed by a factor of cot'^ /3 and, whilst the contribu- 
tion to C(o is enhanced by a factor of tan^ P, it is suppressed by a factor of m^. Including 
BLO effects can alleviate these suppression factors. However, as the Wilson coefficients are 
suppressed by a factor of the overall effect is rather small. 

The LO contributions to Cs,p induced by neutral Higgs penguins and box diagrams 
have been calculated in |2^]. For completeness we present them here 



m 



tan^ f3 y lo^ 



y 



4m^ sin^ Ow 1 - 2/ ' 



<5^-C 



(0)/ 

S,P 



.(0) 
S,P 



(5.9) 



where y = / Tn'jj+ . 

BLO effects can be included by using the couplings presented in subsection |3]^ when 
calculating the matching conditions. The largest correction induced by using these cou- 
plings is attributable to the factor of (1 -|- €3 tan/3) that accompanies the right handed 
coupling of the charged Higgs. This factor typically acts to reduce the charged Higgs con- 
tribution relative to the LO prediction. The GFM corrections to the vertex can act to 
replace the factors of K^^Ktb that characterise flavour change in the MFV contribution 
with the flavour violating insertions ( ^.18D ~( ^T9 ). 

The contributions that arise due to the corrected neutral Higgs vertex proceed via the 
penguin diagram shown in Fig. |l^. Using ( 3.29| ) it is relatively easy to obtain the dominant 
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contribution arising from chargino exchange in the hmit of MFV 



Cs,p 



MFV 



m^rn^At tan^ (3 



H2 {Vun^Vui) 



4m^m^/x sin^ ^h/ (1 + ^3 tan (3) (1 + eg tan 0) 



(5.10) 



denotes the mass of the pseudoscalar Higgs, whilst we decompose the MFV and GFM 
contributions in a similar manner to ( 3.28| ). The most striking aspect of this contribution 
stems from the factor of tan^ /? that appears in the numerator of ( ^.10 ) coupled with a 
relatively weak dependence on the underlying SUSY mass scale. It is therefore possible 
in SUSY models with large tan f3 that large contributions to Cs,p can occur even if the 
sparticle masses are O (1 TeV). (Provided, of course, that the Higgs sector does not decou- 
ple, too.) The BLO contributions in ( |5.1[1| ) are contained in the factors of (1 -l-estan/?) 
and (1 + eg tan/3) that appear in the denominator. In the limit of degenerate sparticle 
masses, for > 0, these corrections tend to reduce the neutral Higgs contribution to Cs^p 
compared with the LO limit ey, — > 0. 

The GFM corrections to the effective neutral Higgs vertex ( 3.30| ) contribute to Cs,p |17| 



nGFM _^ tan^/j 



d-RL TT t ^ „ \xd 



+ 



(l + e^tanp) nib ^ ^ 

d 

(1 + e, tan /?)(! + £3 tan /3) 7^ ' ^'^^ ' ^'^^ ^^^^ 



(5.11) 



and the primed coefficients 



^0 , nGFM _ 4a, uruf, tan^/? 

3a ml KtbKl 

1 

+ 



(1 + estan/?) m-g nis 



(1 + e3tan/3) 



(5.12) 



The contributions arising from the insertions (5^^ and are modified in a similar manner 
to the MFV contribution ( ^.1[)|) and for > undergo the familiar reduction once BLO 
effects are taken into account. It should be noted that once one proceeds beyond the 
approximation of setting electroweak couplings to zero, and uses the substitutions gathered 



in subsection 3.3 an additional contribution, proportional to the insertion, ^^.l ai"ises 



5^ v.. 



EW 



a=l 



tan^/? 



2^/2mwm\ KtbK^^ sin P sin^ 6w "'^^ 



vlH2[yl^,yl^)5i^ (5.13) 



This LO correction tends to interfere destructively with the dominant gluino contribution 
given in ( |5.11| ) and is typically the largest contribution attributable to the insertion df^j^ 
once one proceeds beyond the approximation described in section ^. 

Turning to the insertions and their appearance is a strictly BLO effect and can 
lead to large deviations from LO results where the contributions arising from the insertions 
accidentally cancel as we have shown in |17]. 
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Figure 11: The contributions to Bg fJ-^l-^^ arising from box diagrams mediated by charginos 
and neutralinos. 

5.2 BLO Corrections to SUSY Contributions 

As the gluino does not couple to the leptonic sector the supersymmetric contributions to 
Bg Ai^/f^^ take place via the box diagrams mediated by chargino and neutralino exchange 



shown in Fig. 11. Including BLO effects when calculating these contributions introduces a 



dependence on the bare mass matrix, through the vertices (|2.34D -( |235|) and (|C^)-(|a^ 



Sources of flavour violation can therefore enter through either the chargino or the neutralino 
contributions. However these contributions tend to scale as tan^ /? and, coupled with the 
underlying dependence on l/Mgugy^ rather than 1/m'^, are rather small when compared 
to the effects induced by the neutral Higgs penguins discussed in the previous subsection. 

5.3 Full Calculation 



During our numerical analysis we shall follow the recipe described in subsection O and use 
the expressions gathered in appendix ^ to evaluate the corrections to the various effective 
vertices. We therefore include higher order terms in the MIA as well as any (subdominant) 
SU(2)i X U(l)y breaking effects. Concerning the SM and charged Higgs contributions 
to the decay we use the matching conditions gathered in |Q to evaluate the NLO gluon 
contribution. The contributions that arise from SUSY boxes are given in ||59| , p2| . 

In conclusion, in this section we have discussed how the BLO effects discussed in 
section ^ alter LO contributions to Bg fi^fj.^. In the electroweak sector, these corrections 
typically manifest themselves as factors of either (1 + ea tan /?) or (1 + tan /?) that act to 
reduce both MFV and GFM contributions to the decay. In addition, we have seen that 
new flavour structures, absent at LO, can appear once BLO effects are taken into account, 
leading to potentially large deviations from LO calculations. 

6. Bs - Bs Mixing Beyond the LO 

The final process that we will consider concerns mixing in the Bg meson system. In a 
similar manner to the neutral kaon and B^ systems mixing can occur between the Bg and 
Bg mesons via AF = 2 loop diagrams. In contrast to the neutral kaon and B^i systems, 
however, the mass difference AM^^ between the physical states formed from the two mesons 
has so far remained unobserved. The best bound provided by experiment is currently 

AM^^'P > 14.5ps"\ (6.1) 



In future, the experiments at the Tevatron intend to increase this limit by 20-30% 
whilst even after a year of low luminosity running ATLAS, CMS and LHCb intend to place 
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limits of 30ps , 26ps and 48ps [g^] respectively on AA/^^. Comparing these limits 
with the NLO Standard Model prediction |65, 66 1 

AM|f = (18.0 ± 3.7) ps-\ (6.2) 

it can be seen that the full range of values allowed by the SM can be probed in a relatively 
short time after data taking has commenced at the LHC. 

The effective Hamiltonian that most generally describes Bg — Bg mixing effects is given 
by [|^, H 

We// = [KttKtf) ' (/.) O, (/i) . (6.3) 

i 

In the SM the only non-negligible contribution is proportional to the operator 

qVLL ^ (b"^^PLS^)(bP-f^'PLS^). (6.4) 

However, in the presence of any source of new physics it is possible to induce additional 
contributions to the operators 

Of"" = (b^j^PLsnib^rPRs"), O^"" = {b^PLs''){WPRS^), (6.5) 

qSLL ^ (l^p^s''){b'^PLS^), Oi^^ = {b''a^^PLs''){Wat"'PLsl^), (6.6) 

as well as the parity flipped operators Q^^^ and Q?^^ that can be obtained by substituting 
Pl with Pn in ( |6.4D and (|6.6|). The mass difference AM^^ may then be evaluated by taking 
the matrix element 

AMB^=2\{Bl\'Heff\B% (6.7) 

where {B^\'Hf.ff\B^) is given by 

{m\neff\Bl) =^mlrmBjB. {Kf*Ktf)'Y.P,Cdm), (6-8) 



niB^ denotes the mass of the Bg meson, whilst is given in (5.7). The coefficients Pi 
contain the effects due to RG running between /Xf and as well as the relevant hadronic 
matrix element for the operator in question. Using the MS lattice calculation [39| the 
coefficients Pi have the form 

PY^^ =0.73, Pf^ = - 1.97, P^^ =2.50, Pf^^ = - 1.02, P^^^ = - 1.97. (6.9) 

where we have taken = 4.25 GeV and = rnt{mt). The coefficients P^^^, etc., may 
be obtained by simply exchanging L and R. One interesting aspect of ( |6.9D is that QCD 
effects act to enhance the contributions arising from the scalar operators relative to the 
SM operator C^^^. 

The new physics contributions to neutral meson mixing have been discussed extensively 
in the literature. The NLO charged Higgs contributions to C^^^, for instance, have been 
known for some time now [^] . Turning to the MSSM, the NLO matching conditions have. 
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Figure 12: The tan^ (3 enhanced contribution to arising from charged Higgs and Goldstone 
boson exchange. 



in a similar manner to the decay B Xgj, only been derived for some special cases in the 
MFV limit Most analyses have therefore focussed on using LO matching conditions. 
For example, in ||7^ the LO gluino matching conditions, relevant for mixing in the kaon 
system, were used alongside the NLO anomalous dimension matrix and lattice matrix 
elements to place limits on the insertions governing non-minimal flavour violation between 



down and strange squarks. This analysis has since been extended to the [|7^ and Bs [74] 
meson systems. In the Bg — Bg mixing system, in particular, large contributions to AM^^ 
of up to 120 ps~^ have been found to arise. In the large tan/3 regime another possibility 
for large contributions arises from the inclusion of the effects induced by the neutral Higgs 
penguin. Although strictly speaking an NLO contribution the corrections arising from 
flavour violation mediated by two neutral Higgs penguins have been shown to vary as 
tan^ (5. Such corrections have been analysed in the context of MFV l7|, ^, Q and 
GFM W^, 61, 32, 17] and can similarly induce rather large corrections to AM^^. 



The aim of this section is to discuss how the inclusion of BLO corrections affect the 
electroweak and supersymmetric contributions to the Wilson coefficients associated with 
the operators (|6.4| )-( |6^ . In particular, we shall discuss, in detail, the contributions that 
arise from double Higgs penguins in the GFM scenario and the dominant BLO corrections 
to the existing LO gluino matching conditions. When presenting our analytic expressions 
we shall use the MIA and follow the approximation outlined at the beginning of section |3|. 
We therefore only include the effects induced by gluino and higgsino exchange and ignore 
additional corrections induced by electroweak gaugino exchange and SU(2)i x U(l)y break- 
ing. One may include the additional contributions that appear once one proceeds beyond 
this approximation by using the substitutions gathered in subsection |3.3| . Finally, we shall 



briefly discuss the application of the recipe given in subsection 4.3 to Bg — Bg mixing. 



6.1 BLO Corrections to Electrov^eak Contributions in the MIA 

Charged Higgs exchange leads to contributions to all of the operators given in ( |6.4D ~( |6.6| ) 
(with the exception of Of^^ and its parity flipped counterpart). BLO effects can be easily 
included by using the matching conditions given in [ij] and the corrected vertices given in 
section |3.2| . As an example, the contribution to C|'^ arising from the diagram shown in 
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Figure 13: The contributions to Bs — Bg mixing arising from neutral Higgs penguins. 



Fig. ^ is given by 



5^ a 



LR 



SmfmbUis tan^ /? 
mj^+m'fy (1 + es tan 



m 



w 



2 ' 



(6.10) 



(The other contributions undergo similar corrections.) The factor of (l + estan/?) that 
appears in the denominator acts, once again, to reduce the contributions relative to a 
LO calculation. The GFM corrections included in the term proportional to A32, defined 
in ( 3.17| ), can remove the dependence on the strange quark mass that appears in ( 6.1C| ) 
and lead to potentially large corrections to the MFV result. 

The corrected neutral Higgs vertex can give rise to large corrections to AM^^ in the 
large tan/3 regime via the double Higgs penguin diagrams shown in Fig. |l^ |75|. The 
corresponding matching conditions are 



SLL 



2GWw {KM' 



E 



32 



32 



SO=HO,hO,A° 

E 



m 



Cf 



32 



50 



32 



m 



so 



(6.11) 



(6.12) 



The contribution to the parity flipped operator Of^^ can be obtained via the substitution 
L ^ i? in ( 6.11 ). Using the corrected neutral Higgs vertex ( [j.29 ) in the limit of MFV it is 
easy to obtain the results for Cf^^ [|l^] 



(^^/fO^^iiXMFV ^ GFmlmt 6Hl67r2)^an4/3 



(l + e,tan/3)^(l + e3tan/3)^ 



(6.13) 



and Ck^ 



MFV 



Gprrn^msmi 



el (l67r2)^tan4/3 



^/2^2^2 



^ (1 + estan/3)^ (1 + estan/?)^ 



(6.14) 



Here we decompose the MFV and GFM contributions to the Wilson coefficients in similar 
manner to ( 4.11] ). We have adopted the notation 



sin^ (a 



m 



HO 



m 



hO 



m 



(6.15) 



AO 
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to represent the interference between the scalar and pseudoscalar contributions. At large 
tan/3, ( |6.15 ) becomes (for m\ > m^) 

= ± (6.16) 

and the scalar and pseudoscalar contributions to Cf^^ approximately cancel, whilst the 
contributions to C^"^ interfere constructively. As has been pointed out in [|76| , |l4|, the 
double penguin contribution to acts to reduce the value of AM^^ , bringing it closer 



to the current experimental limit (3.1). The factor of that appears ( |6.14|) suppresses 
the contribution somewhat, however it can still be of the order of 50% of the Standard 
Model contribution in certain regions of parameter space even once the current limits on 
Bg A*"*"/^" 3-^6 taken into account |14]. 

Including the effects of GFM when evaluating the contributions to the Wilson coeffi- 
cients Cf^^ and can lead to far more varied effects. The full expressions for C'l'^ and 
Cf^^ can be obtained by using ( |3.28|) - (^.30|) , however, as they are rather complicated, let 
us highlight the phenomenologically most interesting terms that appear. If the insertion 
d'^i^ is non-zero the contribution to becomes 

X LL TD.Tn ^^ttX% (l67r^)^ tan^ /3 , .2 

^H^'c^Ry = 2 ' 2 -^-^^fe) • (6.17) 

^ 27r2g2 (/^*^Kt,)^(l + e, tan /?)2 (1 + 63 tan /?)2 ^ ^ 

In a similar manner to the MFV correction ( |6.14 ) , the contribution acts to reduce the value 
of AM^^. Once again, however, is suppressed by a factor of the strange quark mass 
that tends to limit the size of the correction to AAf^^. The insertion also contributes 
to via the diagram where one penguin is mediated by gluino exchange and the other 
is mediated by chargino exchange 



)H c^R] = , \ ^ ^ i J^+SU. 6.18 

/ 47r2m2j,K;,i^j,(l + e, tan /3)2 (1 + 63 tan /?) 

In this case the contribution can increase or decrease the value of AMb^ depending on the 
sign of Sfp,. However the factor of that appears in ( |6.18| ) means that these effects are, 
once again, rather small. 

It is possible to avoid factors of the strange quark mass that appear in the above 
expressions for by considering scenarios where the right handed Higgs coupling that 
appears in ( 6.12| ) does not feature such a suppression. This occurs if either <5^^ or (5^^ are 
non-zero. If we consider the diagram where one Higgs penguin is mediated by chargino 
exchange and the other by gluino exchange we have, for non-zero (5^^ [17| 

^/^o^i^xRR = "^b^t eRReyx^^ {iGn^ftan'^P 

2 ) Air'^ml.Kl^Kts (1 + tan /3) (1 + £3 tan (if ''''' 

A similar contribution is possible for the insertion 

{6^ C^^] = ^ \ o^^RT- (6.20) 
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As one penguin is mediated by chargino exchange, ( 5.19| )-( |6.20| ) feature only a linear de- 
pendence on df^p^ and Sfj^. Large positive or negative contributions to AMb^ are therefore 
possible, depending on the sign of the insertions. 

If S'li^ or are non-zero, in addition to either 5^^^ or the diagram involving two 
gluino mediated penguins becomes viable. For example, if S'lp^ and (5^^ are non-zero we 
have the contribution 

^HO^Li?\LL+RR ^ ml ^Rn^LLX-^,^-,^ (l67r^)^an^ (3 ^^^^ 

27r2g2(i^*^K^^)2 (i + e^tan/?)(l + e3tan/3)=^ ''''''''' 



Similar contributions arise for the remaining combinations LL+RL, LR+RR and LR+RL. 
The tan^/3 dependence of (6.21), coupled with its dependence on the strong coupling 



constant present in the factors oi en and eRR, can lead to large corrections to AM^^ in 
the large tan/3 regime. We should briefly mention here how this result compares to that 
presented in |^^. As discussed in |jT^ the authors of |62| omit the effects that arise when 



considering the GFM contributions to the bare CKM matrix (|3.42| ). Once one takes into 
account such contributions, one of the factors of (1 + estan/?) that appears in Eq. (5.8) 
of [^] is replaced by a factor of (1 + e^tan/?). Taking this correction into account, our 
results agree. 

A common feature of all of these corrections lie in the factors of (1 + 63 tan /3) and 
(1 + e^tan/?) that appear in the denominators of all double Higgs penguin contributions. 
These factors represent the resummation of tan/3 enhanced effects and act to reduce the 
contributions for > 0, ^4^ < compared to calculations where resummation is not taken 
into account. 

6.2 BLO Corrections to SUSY Contributions in the MIA 

The super symmetric contributions to Bs — Bs mixing proceed via box diagrams mediated 
by gluino, chargino and neutralino exchange. Unlike the decay B Xs'y however, sizeable 
effects due to these contributions are often limited to regions where the squarks and gluinos 
are relatively light O (500 GeV). 

As Bs — Bs mixing is a AF = 2 process, the gluino and neutralino contributions, at 
LO, feature combinations of two insertions. The LO gluino matching conditions have been 
discussed, in the context of the MIA, in [73|. It is relatively easy to modify them to take 



into account the BLO effects discussed in section ^ by using the recipe discussed at the 
end of subsection [4.2| . The gluino contribution to the operator C'l'^ is given by (to second 
order in the MIA) 



+ 



9a2 {K^^Kts) rni [ (1 + e, tan /3) (1 + e3tan/3) 
42fW (x~) + (^6 + 11^ eLLeRRtan^P \ [2] , , 



mi (1 + estan/3) (1 + e3tan/3)' 



(>LL^RR 



(6.22) 
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The functions f-'"^^ (xj) are given in appendix A. 4 and x^is defined in (4.13). Once again, 



we see that the contributions arising from the insertions dfj^ and 6j^j^ are modified in the 
same manner as the super symmetric contributions to B ^ Xgj (namely ( 4.15| ) and ( [4.16| )) 



that appear in section 4.2. The inclusion of BLO effects therefore tends to reduce the 



overall contribution to arising from these insertions. The BLO term proportional to 
^LL^RR^ on the other hand, tends to interfere with terms that arise at fourth order in the 
MIA and on the whole tends to be rather small. In addition, the combinations (^^^(5^^ 
and (J^/jiJ^/j also appear once BLO effects are taken into account. Due to their tan^ /? 
dependence, they can viably compete with the corresponding terms that arise at higher 
orders in the MIA. Contributions arising from MFV corrections to the bare mass matrix 
can also appear and play the role of Sf^ and insertions. (Contributions where the 
MFV corrections play the role of (5^^^^ and (5^^ insertions are suppressed by factors of m^.) 

Turning to the contributions arising from chargino box diagrams, in MFV the domi- 
nant behaviour at LO arises from contributions to the Wilson coefficient C^^^. At large 
values of tan (3, however, large corrections to Cf 2 are possible that interfere destructively 
with C^^^ and can reduce the contribution that arises from chargino box diagrams to 



AMb^ |14|. The inclusion of BLO effects in MFV, however, tends to reduce this cancel- 
lation somewhat; whilst the Wilson coefficient C^^^ remains virtually unaffected by BLO 
corrections, the BLO corrections to Cf 2 introduce factors of (1 + £3 tan /3) that, in the 
phenomenologically favoured region fi > and At < 0, act to reduce the contribution to 
the Wilson coefficient. These effects therefore reduce the cancellations that occur between 
the two Wilson coefficients, and can lead to an enhancement of the contributions that arise 
from chargino box diagrams [14|. 



In the GFM scenario, the occurrence of the bare mass matrix in the chargino vertex 2.33 
can play the role of one (or both) of the factors of Kts that mediate the flavour change in 
the leading order matching conditions. A dependence on the flavour violating parameters 
can therefore appear once one proceeds beyond the LO. In particular, the factor of 
nis, that features in the matching condition for C|'^ in MFV, can be bypassed in the GFM 
scenario if either 5f^^ or Sf^^ are non-zero. A similar effect occurs for the Wilson coefficients 
Cff^, where non— zero (5^^ or Sf^j^ can lead to additional contributions to the coefficient. 
It is therefore possible to enhance, or decrease, the cancellations that occur between C^^^ 
and the remaining Wilson coefficients depending on the sign of the insertions Sxy- 

The appearance of the bare mass matrix in the neutralino vertex also introduces ad- 
ditional BLO corrections that can modify the LO contributions. In a similar manner to 
the contributions to the decay B Xgj, these corrections typically manifest themselves 
as factors of (1 -|- €3 tan/3) or (1 -|- tan/3) and tend to reduce the contributions compared 
to a LO analysis. 

Let us finally comment on the RG running of the SUSY corrections. The six flavour 
anomalous dimension matrix required to evolve the Wilson coefficients from the SUSY 
matching scale to the electroweak scale was given in [S8|. If we consider the running of the 



-44- 



coefficient C^^ we have 



vr 



25C|'^ {^1SUSY) - ^dC^"^ ifiSUSY) 



(6.23) 



In contrast to the decay B ^s7, the RG evolution of the Wilson coefficients, from ^SUSY 
to //vKi a-cts to increase the coefficients with respect to an analysis where RG running is 
ignored. A similar effect exists in the SLL sector. In the VLL sector, however, RG running 
acts to decrease the Wilson coefficients in a similar manner to 5 ^ Xs^. These effects can 
therefore act to enhance the cancellations between contributions to the VLL sector and the 
remaining operators. 

6.3 Full Calculation 

Let us now discuss the calculation we perform in our numerical analysis. Following the 
method discussed at the end of section ^ we evaluate the SM matching conditions to NLO 
using the matching conditions originally given in p5[ |. For the charged Higgs contribution, 
we use the corrected vertices given in appendix |^ when evaluating the LO matching con- 



ditions given in |14]. When evaluating the SUSY contributions we use the LO matching 
conditions collected in (transformed into our operator basis). They are then subse- 
quently evolved from the SUSY matching scale to the EW scale using the NLO anomalous 
dimension matrix given in [|6^]. 

In summary, in this section we have discussed the effects of all the dominant BLO 
corrections to Bg — Bg mixing. In particular, we have considered all the contributions 
that arise once one takes into account the double Higgs penguin contribution. In the 
region where tan /3 and the sparticle masses are large, for example, these contributions can 
dominate the corrections to AM^^ that arise from new physics, due to the non-decoupling 
property of the corrections to the neutral Higgs vertex. Turning to the supersymmetric 
contributions to the process, we see that, although they play a more minor role compared 
to the neutral Higgs contributions in the large tan (5 regime, the overall effect of BLO 
corrections is to reduce the contributions arising from gluino exchange. We have also seen 
that, in contrast to the decay B — > Xs7, RG evolution between fisusY and fiw can act to 
enhance the contributions of the coefficient in particular. 

7. Numerical Results 

Before proceeding with our numerical results let us first define our parameterisation for the 



soft sector. We treat the soft terms (2.12), defined in the physical SCKM basis, as input. 
For the diagonal elements we set 

("^III),, = ml «^^)^^ = ml 6u, (m|,^)^^ = A, {m,\^ , (7.1) 

The off-diagonal elements are related to the parameters S'Ly via the relations defined in 



(2.18)-(2.19). The soft terms in the up-sector are defined analogously. As inputs for 



-45- 



o 

1—i 

X 




Figure 14: The panel on the left depicts corrections to the Wilson coefficient Cy induced by 
gluino, neutralino and neutral Higgs exchange for MFV. The panel on the right illustrates the 



behaviour of the approximate result (i-12) compared with the full numerical calculation. The soft 
sector is parameterised as follows: rriq = 1 TeV, mg = ^/2mq, At = — 500GeV, ttia = 500 GeV, 
/i = 500 GeV. 



the Higgs sector we take tua, ^ and tan/3 and use FeynHiggs 2.2.8 [|77[ to determine the 
remaining parameters. For the majority of this section we will only vary one (5^y at a 
time unless stated otherwise. Finally the gaugino soft terms M\ and are related to the 
gluino mass m-g via the usual unification relation. 

Let us briefly specify the abbreviations used to denote the various approximations used 
in this section: 

• LO: a calculation that does not feature the resummation procedure described in 
section ^. We do, however, include the LO effects that arise from RG evolution from 
HSUSY to fiw (in contrast with p6[|); 



g-BhO: the approximation used in where only gluino contributions were taken 
into account in the resummation procedure; 

• BLO: the results obtained using the full expressions included in appendix 
7.1 B ^ 

A detailed analysis of the BLO effects relevant to B ^ Xsj was presented in and we 
shall therefore tend to focus on the additional effects induced by the inclusion of electroweak 
contributions. We shall also compare the approximate expressions, gathered in section 
with our complete calculation. 

Before focusing on the GFM scenario, let us briefly consider the effects induced by the 
bare mass matrix and corrected Higgs vertices in the case of MFV. The panel on the left 



of Fig. 14 illustrates the contributions mediated by gluinos, neutralinos and neutral Higgs. 
As is evident from the graph, all three contributions are rather small and, compared with 
the corrections to C-j induced by charged Higgs or chargino exchange, are of the order of a 
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Figure 15: The corrections to tire Wilson coefficient Cy (evaluated at the scale hsusy) induced 
by gluino exchange (on the left) and chargino exchange (on the right). (5^^ = 0.01 and Sf^j = 0.2 
in the top and bottom panels, respectively. The soft sector is described by the same parameters as 
Fig. |lj. As stated at the beginning of section 0, LO is used to denote a calculation that does not 
take into account the BLO effects discussed in section 0. ^-BLO denotes the approximation used 



in 1 15 where only gluino contributions are taken into account in the resummation procedure. BLO 
is used to denote the results obtained using the expressions collected in appendix p[ 



few percent. However, it should be noted that the gluino and neutralino contributions are 
both larger than the neutral Higgs contribution in this particular region of parameter space. 
The panel on the right of Fig. |l^ depicts the approximate expression ( [4. 121) alongside the 
result taken from our numerical analysis. As is evident, ( [4.12| ) describes the behaviour of 
the full numerical result rather well for this choice of parameters and only differs by about 
5% from the result of a full numerical calculation at large tan /?. 



Now let us turn to the GFM scenario. As stated in sections 4.1 and 4.2, it is possible, 
once BLO effects are taken into account, that large cancellations can occur between the 
various supersymmetric contributions. The top two plots in Fig. for example, illustrate 
the cancellations that occur between the chargino and gluino contributions for non-zero 
Sfj^. The dominant BLO chargino contributions to Cy stem from the appearance of the bare 
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Figure 16: Plots comparing of the results of our full numerical calculation with the analytic results, 
derived in the MIA, given in section ^. The gluino contribution is shown in the panel on the left, 
whilst the chargino contribution is shown in the panel on the right. = 0.01 in both figures 

whilst the remaining parameters are described in the caption of Fig. |l^. The MFV contribution to 
the Wilson coefficient in question has been removed in both panels. 



mass matrix in the modified chargino vertex ( |2.34D -( ^.35 ). The inclusion of electroweak 
effects, as indicated by the absence of a term dependent on ey in the terms proportional 
to in ( 4.15| ) and ( [4.20 ), is rather small. The lower two plots in Fig. [l^ show the 
contributions in the primed sector for non-zero Here we see a similar cancellation 

between the BLO effects arising from the chargino contribution and gluino contribution 
to C7. The effect of including electroweak effects is, however, larger than the case of the 
LR insertion due to the appearance of (1 -|- €3 tan/?) (rather than (1 -|- eg tan/3)) in the 
denominators of ( [4.17 ) and ( 4.22| ). However, one can see from the two graphs that the 
increase in the gluino contribution tends to be compensated by a similar correction to the 
chargino contribution. As such the overall effect on the branching ratio is rather small. 

Let us briefly discuss how well the MIA expressions presented in section |^ describe the 
results of our numerical analysis. The two panels in Fig. |l^ show the approximate expres- 
sions for the gluino (4.15) and chargino ( 4.2[1| ) contributions, alongside the results of our 
full numerical analysis, performed in the mass eigenstate formalism where flavour violation 
is communicated via the unitary matrices ( |2.17[ ). As is evident from the graph, in this 
region of parameter space at least, the agreement between the approximate expression and 
that of the full calculation is rather good (within 10%). We have checked that expressions 
for the chargino and gluino contributions that arise from the remaining insertions also tend 
to agree within a similar accuracy. 

Now let us consider the effect of including electroweak corrections when computing 
the corrected charged Higgs vertex (3.14)-(3.15). The panel on the top-left of Fig. 17, 
illustrates the effect of including BLO corrections when calculating the charged Higgs me- 
diated contribution to C7, for non-zero (5^^. If df^ is positive, the GFM corrections to the 
charged Higgs vertex tend to interfere constructively with the MFV contribution, reducing 
the charged Higgs contribution to C7 compared to a LO analysis. On the other hand, if 
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Figure 17: The corrections to ttie Wilson coefficient Cy arising from charged Higgs exchange. In 
the top two panels (5^^ = 0.2 whilst in bottom panels (5^jj = 0.2. LO, BLO and g are defined in 
the caption of Fig. |l^. In the panels to the right our analytic results are compared with those of 
our complete numerical analysis. The MFV contributions to Cj and arising from charged Higgs 
exchange have been removed in the panels to the right. The abbreviation "Approx." denotes the 



results gathered in section 4.1, whilst "Impr. App." is used to denote the same results supplemented 



by the substitutions gathered in subsection B.S . The soft sector is described by the same parameters 
as Fig. |l^. 



is negative, the MFV and GFM corrections to the LO contribution interfere destructively 
and can therefore lead to an enhancement of the charged Higgs contribution to the decay. 
As Ai, is set equal to zero, the bulk of the difference induced by the inclusion of EW contri- 
butions in the figure results from the correction to the right handed coupling ( p.l5 ) and the 
corrections that arise from the additional gaugino mediated electroweak effects discussed in 
subsection 3.3. This is confirmed in the top-right panel that illustrates the behaviour of the 



approximate expression ( [4. 61 ) compared with the results of our full numerical calculation. 
As is evident from the figure, the agreement between the curves is rather poor and the 
approximate result where the electroweak couplings gi and §2 are ignored provides a 30% 
overestimate of the beyond leading order effects. The origin of this discrepancy stems from 
graphs featuring gaugino exchange. This result is not specific to the GFM scenario and 
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such effects liave been discussed before, in the context of MFV in |14]. Once one takes into 
account these effects, by applying the substitutions found in subsection 3^, the agreement 
between the numerical and approximate results improves dramatically. As is evident from 
the line depicting the improved approximation in the top-right panel of Fig. 17. 

The lower two panels in Fig. depict the contributions to the primed coefficients 
arising from BLO corrections to charged Higgs exchange for = 0.2. In a similar 
manner to the insertion S'l^, the majority of the difference between the effects considered 
in (the line labeled ^^^BLO) and the complete calculation presented in this analysis, 
stems from the destructive interference between the e<j and ey terms that appear in the 
denominator of (3.15) and the additional electroweak corrections to the (3,3) element of 
the left-handed charged Higgs vertex. Both of these effects act to increase the charged 
Higgs contribution attributable to RR insertions. This situation is also evident in the 
lower right panel, where the inclusion of the additional electroweak corrections described 
in subsection |3.3| roughly doubles the accuracy of the approximate expression. 



The contributions due to new physics for each insertion are shown in Fig. 18. As is 
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Figure 19: The variation of BR(B — 
sector is pararaeterised as follows rriq 



Xs7) with the flavour violating parameters The soft 



ITeV, 



^500 GeV, ruA = 500 GeV, 



Ai = 500 GeV and tan/3 = 
shown in all four panels. 



50. A broad region in agreement with the current experimental limit is 



evident from the plots for the LR, RL and RR insertions, the effect of BLO corrections 
tends to be rather large. For example, in the case of the LR insertion (the top-right 
panel in Fig. 18), the reduction of the gluino contribution with increasing tan /?, coupled 
with its destructive interference with the chargino contribution, can dramatically alter the 
behaviour of the Wilson coefficient for even moderate tan /3. Similar effects occur for the 
RL and RR insertions. For the LL insertion the difference between the LO and BLO 



calculations shown in the top-left panel of Fig. 18 is rather slight. This is because a tan /? 
enhanced correction, proportional to 5^^;^, appears at LO in the chargino Wilson coefficients. 
It is therefore rather difficult for the BLO corrections to the chargino contribution to play 
as large a role as they do for the other three insertions. 

Finally, let us consider the combined effect these contributions have on the branching 
ratio. Fig. ^ illustrates the variation of the branching ratio with the flavour violating 
parameters As is evident from the figure, BLO effects can be significant for all four 

insertions. Contributions due to flavour violation in the LR and RL sectors, in particular. 
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Figure 20: The variation of BR(i3 Xs^) witli the tan/3. The soft sector is parameterised as 
follows niq = ITeV, = VSm^, A = -500 GeV, ruA = 500 GeV and /i = 500 GeV. A broad 
region in agreement with the current experimental limit is also shown when applicable. 



can undergo reductions by up to a factor of two compared to a LO analysis. Turning to the 
inclusion of electroweak effects, we see that although such corrections can effect the Wilson 
coefficients by up to 30%, the overall difference between the approximation used in |jl6[ 
and the full calculation used in this analysis is rather small. This is primarily because the 
origin of the large discrepancy, between the LO and BLO calculations, is mainly due to the 
cancellation between the gluino and chargino contributions, that arise for each insertion. As 
both the chargino and gluino contributions undergo similar corrections, once electroweak 
effects are taken into account, the overall effect on the branching ratio tends to be rather 
minor. 



Fig. 20 depicts the variation of the branching ratio with tan (i. All four panels exhibit 
the focusing effect first described in |l^ where BLO corrections act to reduce the LO 
result such that the SM result ( [4.2D is preferred. Once again we see that for the insertion 
(5^^ BLO corrections are rather small due to the presence of the LO correction to the 
chargino contribution in ( 4.18| ). For the remaining insertions we see that the BLO effects 
described in this paper can significantly alter LO corrections for tan/3 as low as 20. 
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In summary, we have seen that whilst electroweak corrections can effect individual 
Wilson coefficients by up to 20%, the overall effect of such contributions is rather small. 



In particular, the focusing effect described in [15, 16| remains even once one includes all 
the electroweak corrections described in this paper. This is unsurprising as the focusing 
effect mainly arises from the combined effect of two contributions. The first arises from 
renormalization group evolution of the Wilson coefficients evaluated at the SUSY scale 
to the electroweak scale, which is naturally independent of electroweak corrections. The 
second is due to cancellations between the gluino and chargino contributions to the de- 
cay. We have seen however that electroweak corrections typically alter each contribution 
in a similar manner, and an increase in the gluino contribution, attributable to the de- 
structive interference of the gluino and chargino contributions to the bare mass matrix, 
is often accompanied by a similar increase in the chargino contribution. Finally, we have 
seen that the approximate expressions gathered in section ^ tend to describe the overall 
behaviour of the super symmetric contributions rather well. The contributions arising from 
the charged Higgs exchange however, often have to be modified according to the improved 



approximation described in subsection 3.3 to obtain a sufficient level of accuracy. 
7.2 Bs ^ fi+fi- 

As discussed in section ^, large corrections to the decay Bg — > lJ-^lJ-~ are possible in the 
large tan P regime due to the contributions of scalar and pseudoscalar operators that are 
proportional to tan'^ (3. Fig. 21 depicts the dependence of our analytic and numerical 
results for the Wilson coefficients Cp and C'p on tan/? for various choices of flavour vi- 
olating parameters. The top-left panel in the flgure shows the contribution due to the 
insertion S'li^. Let us point out that when one only uses expression (5.11) to calculate 



the contribution to Cp, one obtains a value roughly 20% larger than the result obtained 
from a full numerical analysis. Once one takes into account the additional LO correction 
induced by gaugino-higgsino mixing ( 5.13| ) (denoted by Impr. App. in Fig. 21) the agree- 



ment improves to roughly 10%. The remaining sources of discrepancy are mainly due to 
wavefunction corrections to the bare mass matrix, which can be as large as 10%, and the 
inevitable limitations associated with the MIA. The contributions due to the insertions 
(5^^ and (5^^ shown in the top-right and bottom-left figures respectively are absent at 
LO as the effects of the insertions cancel [0]. However, once BLO effects are taken into 
account, a dependence on the insertions is reintroduced due to their appearance in the 
bare mass matrix (|3.8| ). (The dependence of (|3.S|) becomes apparent once one recalls 
that (5^/2)23 ~ ('^^i)32') Comparing the approximate and numerical results for the BLO 
corrections we see that they typically agree with one another very well. (Unlike the df^^^ 
insertion, contributions due to wavefunction corrections are typically rather small as they 
appear at second order in the MIA and tend to be suppressed by factors of m^.) Finally, 
the bottom-right panel shows the contributions that arise for non-zero 6j^^ where, once 
again, the analytic and numerical results agree with one another rather well. 

With these results in mind let us now consider the overall effect of such corrections on 
the branching ratio of the decay. Fig. ^ depicts the dependence the branching ratio on 
the various sources of fiavour violation in the squark sector. The two graphs depicting the 
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Figure 21: The tan/3 dependence of the Wilson coefficients Cp and C'p. In the top two panels 
0.2 and Sfj^ = 0.01, respectively, whilst in the lower two (5^^ = 0.01 and (5^^ = 0.2. In all 



four panels the MFV contributions to Cp have been removed. The soft sector is parameterised 
as follows rriq = 1 TeV, rrig = V^niq, An = -500 GeV, toa = 500 GeV, = 500 GeV. The 
abbreviation "Approx." is used to denote the contributions arising from formulae (|ti])-(|t|) 
whilst "Imp. Appr." is used to denote the calculation that includes the additional electroweak 



contribution (5.13) 



variation with and (^^^ show the characteristic suppression associated with the BLO 
factors of (1 -|- €3 tan j3) and (1 -|- tan 0) that appear in the denominators of the Wilson 
coefficients This suppression can loosen the bounds placed on these insertions 

by the Bg fJ-^fJ-~ constraint. The panel depicting the variation with (^^^ displays a 
larger dependence on the the gluino mass than the panel featuring the insertion 6j^j^ as 
the gluino mass not only features in the gluino contribution, but also in the corrections 
that arise once one includes gaugino-higgsino mixing. (We remind the reader that we 
assume that the mass of the wino and the gluino are related.) Let us briefly comment 
however, that in contrast to B ^ ^s7) the differences between the BLO and LO results 
for MFV and GFM calculations tend to be rather similar. This is because, now, the only 
dominant contribution to the decay is via the neutral Higgs penguin and BLO effects 
tend to be limited to the factors of (1 -|- tan /3) and (1 -|- €3 tan /?) that accompany the LO 
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Figure 22: The dependence of BR{Bs M^M") on Sj^y for = 1 TeV, Au = -500 GeV, 
TTiA = 500 GeV, /i = 500 GeV and tan /3 = 50. The Ught blue (hght grey) and red (dark grey) bands 
depict the effect of varying rrig between mq/\/2 and y/2m,q. The pubhshed D0 and prehminary 
CDF hmits for the decay are shown. 



matching conditions. Turning to the and 6"^^ insertions, large deviations from a purely 
leading order calculation are possible due to the reappearance of the insertion in the Wilson 
coefficients (|5ll|)-(|5l^) once one proceeds beyond the LO Finally, the asymmetric 
nature of the d'^^^ and curves arises as these contributions interfere directly with the 
MFV contribution and it is therefore possible to induce quite large cancellations. For dj^j^ 
and on the other hand, direct interference with the MFV contributions is generally 
not possible. Cancellations with the Wilson coefficient Ca when calculating the branching 
ratio ( p.5| ) can occur, however, and lead the minima of the curves to deviate slightly from 
MFV. In general however, in a similar manner to i? ^ Xs'j, the overall effect of these 
insertions is to increase the branching ratio with respect to the MFV result, independent 
of the sign of the insertion. 

The dependence of the branching ratio on tan (3 and the pseudoscalar mass tha is shown 
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Figure 23: The dependence of BR(i3s jjL^ (x~) on tan /3 (the top two panels) and (the bottom 
two panels). The soft sector is parameterised as follows: mq = ^/2mg = 1 TeV, A„ = —500 GeV, 
H = 500 GeV. In the top two panels tjia = 500 GeV, in the lower two panels tan/3 = 40. The 
published D0 and the preliminary CDF bounds are shown when appropriate. Only one LO curve 
is shown in the panels on the right as the dependence on the insertion vanishes at LO. 



in Fig. 23 . As discussed in section g the scalar and pseudoscalar contributions to the decay 
can lead the branching ratio to vary as tan^ /3. As is evident from the top two panels, values 
of BR(fis — > /i+zi") approaching (or even exceeding) the current experimental limit are 
possible, even for TeV scale sparticle masses, if tan /? is large ~ 40. A strong dependence 
on niA is also apparent in the lower two panels. Both figures illustrate the reductions 
associated with BLO calculations for the LL insertion and the new effects that appear 
beyond the leading order for the RL insertion. 

In summary, we have seen that the approximate formulae gathered in section || seem 
to describe the results of our numerical analysis rather well (within 10%) especially once 
one considers the approximations involved in their derivation. For the LL insertion the 
additional electroweak corrections, described in subsection W^, typically act to reduce the 
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Figure 24: Contributions to AM^^ vs. tan j3. The combined effect of box diagrams mediated 
by SUSY particles is shown in the left panel whilst the right panel shows the double penguin 
contribution. The soft sector is parameterised as follows: = 500 GeV, At = —rriq, nig = ^/2mq, 
ji — ra-q and toa — 300 GeV. The flavour violating parameters are — (5^^ = 0.1. 



correction to Bg — > ^^/^^ by up to 20% compared with calculations performed in the limit 
of vanishing electroweak couplings. This reduction coupled with the resummation of tan (3 
enhanced effects can relax the contribution due to 5'^^ by roughly 60% compared with a 
naive LO analysis in which only the effects of the gluino contribution to the neutral Higgs 
vertex are taken into account. 

7.3 Bs - Bs Mixing 

Turning now to Bg — Bg mixing, here we shall not compare the expressions gathered in 
section ^ to those of our numerical analysis as our approximations for the effective Higgs 
vertex have been discussed in the previous subsection. 

As discussed in section contributions to AM^^ in the GFM scenario stem from box 
diagrams mediated by SUSY particles and charged Higgs and, in the large tan/3 regime, 
from double penguin diagrams. The panel on the left of Fig. ^ shows the tan (3 dependence 
of the contributions to AM^^ arising from box diagrams mediated by SUSY particles. The 
main difference between the two curves at low tan/3, originates from the use of the NLO 
anomalous dimension matrix to run the BLO calculation from the SUSY matching scale 
to the electroweak matching scale. At large tan/3 the interference between the dominant 
gluino contribution and the BLO corrections to the chargino and neutralino contributions 
acts to reduce the overall contribution to AM^^ further. As is apparent from plot however, 
the overall correction tends to be only of the order of five to ten percent. 



The double penguin contributions are depicted on the panel on the right of Fig. 
Here we see a rather more dramatic difference between LO and BLO calculations and at 
large tan (3 it is possible that BLO corrections can lead to a reduction of LO effects by up 
to a factor of three. 

Let us briefly comment on which values of tan /3 should analyses, that typically only 
feature LO matching conditions for the gluino contributions to Bg — Bg mixing, include 
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Figure 25: The dependence of AAIb^ on Sj^y foi' the same parameters as Fig 



the double Higgs penguin contribution. As is evident from Fig. the double penguin 
contribution completely dominates the behaviour of AM at large tan (3 and can become 
as important as the LO result for tan/5 as low as 20. We have checked that a similar 
situation arises for non-zero LR and RL insertions. 

With these results in mind let us consider the combined effect of all the beyond Stan- 
dard Model corrections in the large tan /3 regime. The dependence of AMb^ on each of the 
flavour violating parameters is illustrated in Fig. |2^. In the top left panel corresponding to 
the insertion 6f^^ one can see a largely quadratic dependence on the insertion, in agreement 
with the analytic result ( 6.17| ). The graph is not centred on = as, in a similar manner 
to the Bg — At"*"/^" graph in Fig. the MFV and GFM contributions to the neutral Higgs 
vertex approximately cancel for Sf^^ ~ 0.08. Turning to the top right panel, depicting the 
dependence of AMb^ on we can see that the overall effect on AMb^ is rather slight, 
for small the effect is mainly linear as only contributes to the left-handed vertex 
that appears in ( 6.121) . It is therefore necessary for the other penguin to be mediated by 
chargino exchange (|6.18|) . For larger values of 5^^, contributions to Cf^^, as well as SUSY 
box diagrams, lead to a quadratic dependence on the insertion to emerge, however, once 
again, the corrections are rather small. 
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Figure 26: Scatter plots showing the correlation between AMb^ and BR(i?s fi^fi^). In the 
panel to the left S'lj^ is varied over the range [—0.8, 0.8] whilst in the panel to the right the insertion 
^RR varied over the same range of values. In both plots niq is varied over the range [500, 1500] GeV. 
The remaining parameters describing the SUSY sector are as follows: /i = 500 GeV, ttia = 500 GeV, 
Au = -500 GeV and mg = 1 TeV, tan^ = 40. The published D0 and the preliminary CDF bounds 
for the decay Bg fJ'^ 1^^ are shown, as well as the lower limit for AMb^- Points compatible with 
B Xgj are highlighted by light blue (light grey) squares. 

The bottom two panels depict the larger effects induced in the (5^^ and (5^^ sectors. 
The linear dependence of the contributions is due once again to one Higgs penguin being 
mediated by chargino exchange and the other by gluino exchange. The only alteration to 
this behaviour arises for very large (5^^ (~ 0.4) where the gluino mediated contributions 
to the left-handed Higgs coupling, that are suppressed by a factor of rus, can become 
important and interfere with the chargino contribution. Finally, let us once again point 
out the large differences between the LO and BLO calculations featured in all four plots. 
For dfj^ and (5^^ we see the characteristic suppression of LO effects that arise from the 
factors of (l + e3tan/3) and (l + e^tan/?) that appear in (|6l^ )-(|6^). These typically 
lead to reductions proportional to factors of two or three, if > 0. In a similar manner 
to the decay Bs —>■ fJ'^^J'', a dependence on the insertions 6'[^ and 5^^, which is absent at 
LO, reappears when BLO corrections are taken into account [|l7| ]. 

Before ending this section let us briefly discuss the correlation between Bg — > fJ-^fi^ 
and AMb^ at large tan/3. As was pointed out earlier in this subsection, the double Higgs 
penguin tends to completely dominate the contributions that arise from new physics in the 
large tan /? regime. It is therefore natural to expect a degree of correlation with the decay 
Bg fJ-^fJ-" , that is also dependent on the neutral Higgs penguin when tan/3 is large. Such 
a situation is illustrated by the scatter plots shown in Fig. Here we see in both panels 
that a large correlation exists between each process. The panel on the left (where only the 
LL insertion is varied) features only one branch as the double Higgs penguin in this case 
tends to only lead to a reduction in AMb^ ■ Varying the RR insertion can lead to reductions 
or enhancements of AM^^ thanks to the linear dependence on the insertion exhibited in 
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Figure 27: The ratio K^^, /K^f vs. (5^^ (the panel on the left) and (5^^ (the panel on the right). 
The soft sector is parameterised as follows: tan (3 — 40, — 1 TeV, ji — 500 GeV, mA = 500 GeV 
and At = —1 TeV. Solid and dashed lines denote nig — \/2mg and m.g = rriq /\/2 respectively. 



the matching condition given in ( |6.19| ). This effect leads to two distinct branches being 
visible in the right panel of Fig. 

7.4 Radiative Corrections to the CKM Matrix 



It was mentioned briefly in subsection 3.4 that the rotation from the bare to the physical 
SCKM basis can induce large contributions to the bare CKM matrix K^'^^ |2^. One 
particularly interesting consequence of this is that the CKM matrix elements K^s and K^t 



could be generated radiatively via GFM effects |16|. Such a situation is illustrated in 



Fig. 27. From both plots it is evident that (and similarly K^^^) can receive significant 



^cb 

corrections due to the presence of GFM in the squark sector. In the LR sector in particular, 

radiative generation of the entire matrix element can typically occur for dfj^ 0.01. 

One could imagine the scenario where corrections to the remaining matrix elements might 
induce similar effects and lead the CKM matrix to be fully diagonal before SUSY threshold 
corrections are taken into account. 



From the top-right panel in Fig. 22 it is also apparent that, due to a cancellation 
between the chargino and gluino corrections to neutral Higgs vertex, the branching ratio 
of the decay Bs fJ-^fJ-^ ^-Iso tends to approach SM like values for negative 5^^- The 
minima of the two curves however tend not to coincide unless eyY^^ ~ — (this typically 
only occurs if At is rather large and negative). 



8. Summary 

We have presented here the first complete analysis that includes the resummation of all 
tan P enhanced BLO effects in SUSY with GFM that contribute to the processes B Xgj, 
Bg Bs — Bg mixing. We have derived analytic expressions applicable in general 

and in the MIA. As such, they match the precision of similar calculations performed in the 
context of MFV. We have provided a recipe for including BLO effects into LO expressions. 
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We have found that BLO effects in GFM can be large. In particular, we have more fully 



analysed the focusing effect that was initially pointed out in 16 1 in the case of 5 — > 



and next shown to also exist in the other two processes |17]. In the phenomenologically 
interesting case of large tan (3 and /i > the focusing effect often leads to a significant 
relaxation of experimental bounds on the soft mass mixings. 

Finally, we have examined radiative corrections to the CKM entries. These can be 
large, or even dominant, due to large LR mixings. We have pointed out a possible, although 
apparently accidental, correlation with the processes analysed here. 

The method presented here is rather general and can readily be extended to include 
CP violating contributions, here BLO corrections can carry additional phases and lead to 
potentially large deviations from LO calculations. One might also consider the process 
B Xsl~^l~ , where the BLO corrections to the Wilson coefficients Cj and Cg (and their 
primed counterparts) may also play a large role. Another possible application would be to 
include the Bd — B^ mixing system where, at large tan /?, one can expect strict bounds the 
insertions (<J^/j)]^3 and {^111)1^ due to contributions to AMb^ unsuppressed by m^. 

The formalism presented here can be applied to data from present and future B- 
factories and hadron colliders in constraining mass insertions and, eventually, in extracting 
information on an emerging pattern of flavour violation in the squark sector ||7q]. 
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A. Loop Functions 
A.l The Functions Hi 

The loop functions Hi that appear throughout the text are given by 

rj r N Xl log Xl X2 log X2 , ^ 

=(l-x,)(x,-X2) + (1-X2)(X2-X0' ^^-'^ 
/73(x„X2,X3) (A.2) 

Hi(xx,X2,X'i,XA) = . (A.3) 

X3 - 



- 61 - 



In the limit of degenerate arguments the functions become 

^2(1,1) = -^, (1,1,1) = ^, if4(l,l,l,l) = -^. (A.4) 

A.2 B Xsj 

The loop functions g* (x) that appear in the charged Higgs matching conditions ( |4.5| )- 
( |47| ) are 

^(2) , ^ X (3 — 5x) X (3x — 2) , 

(x) = — ^ + ^ ^ log X, 

(2) _ x(3-x) ^ 1 

The functions i?]"^'^' {x) that appear in section |4.2| are given by the following expres- 



sions. The chargino contribution features the functions 

„n] , , -Sx^ + 2x , -Sx^ - 5x + 7 

iil (x) = IT log X H ^ — , 

^ 6(l-x)^ 36(l-x)3 

jj[8] .X X -x2 + 5x + 2 

//{ (x) = J logx H o— , 

^ ^ ^ 2(l-x)^ 12(l-x)^ 



^ f7l , s -3x2 _^ 2x ^ -5x2 _^ 3^ 

-Ha (x) = 0- logx H ^, 

2 3(l-x)^ 6(l-x)2 

-x^ + 3x 



(x) = — — ^logx + 



X 



(1 — x) 2 (1 — x) 

The functions relevant to the neutralino contribution are given by 

hI\x) = -\Hf {x) , Hf{x) = Hf{x), 

i/f (X) = -\ [nf (X) + i) , Bf (X) = Bf (x) + \. 

The functions relevant to the gluino contribution are 

The loop functions /] ' (x) and J\ ' (x) that appear at higher orders in the MIA are 
related to the functions B\ ' (x) via the relations 
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A.3 Bs fi^fi' 

The function that appears in (|5.8|) is given by 



, , , X , (x3 - 6x2 + 3x + 2) 

fz{x) = — -4 logx + — -4 . 

4(x-l) 24(x-l) 

A.4 Bs - Bs 

The functions that appear the gluino matching condition ( |6.22 ) are given by 
^[11,, x2(x + 3)^ (-17x3 + 9x2 + 9x-l) 



^ ' ' (x-1)^ 6(x-l)^ 

^[21 , , 2x2 (x + 1) ^ (-x3 - 9x2 + + l)^; 

/- (2;) = K- logx H E . 

f ^ ^ (x-1)^ 3{x-lf 

A. 5 Passarino Veltman Functions 

The Passarino- Veltman functions that appear in appendix ^ have the following form 
Bo (2;, y) = - log ^ + IH — log -, 

1 1x1 X y^ — 2xy y 
5i X, y) = --1] + - log ^ - - - — + — -2 log -, 

2 2 /i^ 4 2 (x - y) 2 (x - y)"' x 

C^o (a;, y, = 7 77 r log - + log -, 

{x-y){y-z) X [x-z)[z-y) X 

^ , n1/^3,x y2 y ^2 ^ z\ 

Coo (x, y, z) = - r/ + - - log — + log - + log - , 

4 \ 2 [X — y)[y — z) x [x — z)[z — y) xj 

where = f + log An — ^e- 

B. Alternative forms for 6'^^^ and 

As discussed in section throughout this analysis we assume that the trilinear SUSY 
breaking terms are not proportional to the appropriate Yukawa coupling. To illustrate how 
our results are altered if we do make this assumption, let us consider a specific example 
with relevance to models with SUSY breaking mediated by either supergravity or gauge 
interactions [^] 

ml^^ = A\^mf^ + m'':P^A\^. (B.l) 
The insertions (5^^ and are therefore equal to 



23 



23 



(B.2) 
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It should be noted that, in addition to the off-diagonal elements of the complex matrices 
A^'^, there is a contribution proportional to the appropriate off-diagonal element of rn"^\ 
If we use the choice ( [B.lD , one can redefine the insertions and 5'^i^ ( |B.2[) to be inde- 
pendent of the bare mass matrix in the following way (ignoring terms suppressed by the 
strange quark mass) 



^III) 22 {^Irr) 33 ^ [^III) 22 [^IrR^ 33 

This definition is largely independent of BLO corrections and would correspond to a input 



value appropriate for the iterative procedure described in section 2.4. The effects due to 
the bare mass matrix then appear in two ways. The diagonal elements result in a factor of 
1/(1 + 63 tan P) that accompanies each factor of (5^^ and 5^^^. The off-diagonal elements, 
on the other hand, may be included in a similar manner to the off-diagonal elements of 
the F-terms that appear in the squark mass matrix by altering with a small (cot /? 
suppressed) correction. The only exception to this rule concerns the corrected charged 
Higgs vertex, where one must perform the substitution 



;;d _^ , \ "^7 22 V / 23 V 7 33 V / 23 
'^LR '^LR I ^ = 

TflP' 1 f Tfi^ 1 
'^'■^^/22 V d,RR)^^ 

in ( |3.21| ). In the limit of MFV [Adi = ^O) ^di? = 0) we reproduce the result presented 
in 0. 

C. Supersymmetric Vertices 

Now let us present a complete list of supersymmetric vertices required for our calculation. 
Throughout this section /, J = 1, . . . , 6, a, 6 = 1, 2, a, /3 = 1, . . . , 4 and finally j = 1, 2, 3. 
The coupling of the gluino to down quarks and squarks is given by 

Cg = d\ {g) [{GdL)j, Pl + {GdR)j, Pr] {d\ , (C.l) 

where {GdL)ji and {GdR)ji are given by 

{GdL)j, = -V2gs {TdL)ji , iGdR)j, = V2gs {TdR)j, (C.2) 

The couplings to up quarks and squarks may be obtained via the simple substitution d ^ u. 



The chargino coupling to down quarks and up squarks was discussed in section 2.3 
whilst the coupling to up quarks and down squarks is given by 

iCuL)aJi = - 92Ual (TdLK^ + -=-^ -C/„2 (VdRm^^^K^) (C.3) 

\ / Jt y/zmw COS p ^ ^ 

{CuR)aJt =^^^^—-Va2 (r,^i^t^(0)t) (C.4) 
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The neutralino couplings between down quarks and squarks are 

^x' = d\{x'')J{NaL)o..HPL + {NdR)^j,PR]{d)^, (C.5) 
where {NdL)jia and {Ndnjji^ are given by 

V2 V-J / V 2r?T.vi/ cos p ^ 

(C.6) 

(^dR)aJi = o ^"1 (rdR)j, - -7= -Na3 i^dLm'^ j . (C.7) 

The matrix diagonahses the neutrahno mass matrix in the usual manner 

N*M^oN'^ = diag (^m^o,. . .,m^(^ . (C.8) 
The couplings to up quarks and squarks are 



(C.9) 

{Nur)^j, = "''""T""" iv»i (r„fl)^, - ^ . ( r„Lmrn ^. . (C.IO) 



2\/25r2 tan 6*14/ ^ 52 „„ {o)t 



Let us now consider the couplings of the W boson that appear in ( D.3 ). The coupling 
to squarks is given by 

W.j = -^(rnLi^rtJ^^. (c.ii) 

The couplings to neutralinos and charginos are given by 

{W^L)aa=92{Nl^Val-^K^Va2^ (W^^) = (^N^^Ual + ^N^zUa^ (C.12) 

The couplings of the Z boson to up and down squarks are given by 

(t. T.t 4 2, 



The couplings to neutralinos and charginos on the other hand are 

[UaiU^i + (cos2 Ow - sin2 Ow) 6ab\ , (C.15) 

VC^ + (cos2 Ow - sin^ Ow) ^ab] , (C.16) 

{K,Np,-K,Nps)- (C.18) 







92 r 


2 cos 9w 




R)ab = 


92 r 


2 cos 6w 






92 . 


L)al3 ~ 


2 COS Ow 




R)a/3 ~ 


92 

2 cos 9w 
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The coupling of the charged Higgs boson to squarks is given by 



) IJ V2sin6'H/ 



IJ 



+ 
+ 



92 



92VS+ 



+ 



+ 



cos/3 v^ """""'^ ^^//j ' sin/3 

\/2cos/3sin/3 V " <^ <iR) ij 

TuLKmd,LR^\j^ + 11* cot/3 {VuRm^^^KT^^ 

{ruRm^,RLKV\^) + /xtan^ (r^^Kmf tpt ^) 



IJ 



+ y(2) 



\/2mvKCOs/3 - 
'■■^^ \/2mw sin /3 



/J 



7JJ 

(C.19) 



where rjg+ = 1,0 The couphng of the charged Higgs to the chargino and neutrahnos is 
given by 



( 



C: 



s+ 

xL 



92 S+ 

y{2) 



cos Ow 



1 

71 



Ua2 (Ki sin Ow + N*2 cos ew) - UaiN*3 COS ew 



r^^Jaa" COS^^y^W 



Ki2 (-/Val sin 9w + -/Va2 COS 9w) + 14i COS 6*1^ 



(C.20) 



. (C.21) 



The couphngs of the neutral Higgs bosons to squarks are 

/c-S°^ iy^^ 2n ( s'> ■ o s"\ f-, , 3-4sin^ g^y . 

\ d Owmw [cos Px(^^^ -sm(3x^^2)) [^+ o.i„2Q__. ^dL^'^]^ 



92 50 
^(1) 



mw cos (5 
92 

2mw cos P 



2 sin"' 9w 



IJ 



(C.22) 



3-8sin2 Ow 



='?5o| - ^ tan^ 6lw'mv^ (^cos /Jxg^"-, - sin/Jxg,) + ^^.^^ 



IJ 



92 ^S" 

mw sin P 



/J. 



92 



2mvi/ sin /3 

where 7750 = 1,1,0,0 and = sin a, cos q, i cos /?, z sin f3 

The couplings of the bosons to charginos and neutrahnos are 



(C.23) 
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i^X-R)ab 



92 ( S°*TT T^* I 5^*77 jr*\ 



(C.24) 
(C.25) 



92 



2 cos 6*^ 

5°* 



6"* AT* 

^(1) ^^a3 



+ ( xfy*A^^3 - xf2;A^^4j (Ki sin^H/ - K2 COS 9w) 



X{2)*K4) (^^1 sin - N*p2 cos 



92 



'P 2 cos 6*^^ 



X^2°)^a4 ) (A''/31 sin Ow - Np2 COS 6*1^) 



+ - 2:^2)^/34) (^^al sin 6*14/ - iVa2 COS Ow) 



(C.26) 



(C.27) 



D. Corrected Vertices 

Let us now present the full analytic expressions used in our numerical analysis. A number of 
the corrections given in this section are divergent and should be renormalized appropriately. 
However, we have checked that numerically these terms are negligible compared to the 
tan (3 enhanced corrections discussed in ^ and tend to play only a minor role. Throughout 
this section i,j = 1,2,3, I,J = 1,...,6, a,b = 1,2, a, /3 = 1,...,4, 5"+ = H^,G'^ and 
50 = H^,h^,A^,G^. Repeated indices that appear in the expressions below should be 



summed over. We have checked are results with those that appear in [14| and find that, 
once one takes into account the different form of Passarino-Veltman functions we adopt 
our results agree. 



The self-energies that appear in the corrected vertices presented in section |2.2| as well 
as the expression for 5ma are given by 



Ed 
mL 



1 



m-g {G*^R)n {GdL)ij C2 (3) (m|, m| 



167r2 



(D.l) 



11 



167r^ 



+ iCdL)aii iCdL)ai,Bi[ml.,mi 



(D.2) 



The various supersymmetric couplings that appear in the above expressions are given in 
appendix whilst the Passarino-Veltman functions are defined in appendix A. 5. can 
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be obtained by substituting L with R in (|D.2| ) the up quark self-energy corrections can be 
obtained by substituting u ^ d. 

The correction AC^ that appears in ( |2.24| ) is given by 



(^^^ =16^1^^^ ^^'"^^^^ H'<'< 

+ 2 {KL)aU {NdL)^jj Coo (m^o ,m|^,mi^) 

+ "^X^a (^xR)aa iKL)aR iCdL)aIj Co (mjo , _ , m| 



"^Xna (^xL)aa iC*L)aM {^dL)ajj Co (^m^o , , 



Coo ("i^S, 
Coo 



(D.3) 



can be obtained by the simple substitution L ^ R. The correction to the left handed 
coupling of the Z-boson ( |2.27 ) is given by 

=li{'^^ ^^S)u iG*,L)n iG,L)j, Coo {ml,ml,ml 
+ 2 {Zs)u iNh)aH iNdL)^j, Coo (m^o ,m|,m|^ 
+ 2 {Z^),j {C2L)aH iCdL)aj, Coo (m^. , , mi^ 



(^x-i)afe ('^rfi)a/j Co (m^ _ , _ , m| J 



"^x°"^x^ (^x°^^)a/3 (^rfL)/3J* (^'ii)aJi ^0 (m^o,m^o ,m~^ 



+ 2 Z 



X" -RJab 



{CdL)bIi {CdL)a 



Act 



+ 2 (^X°i?)a/3 (^rfL)/3J* (^di)aJi CqO ("^x° ' ""x)^ ' "'IJ " i 

AC^ may be obtained in a similar manner to AC^ . 

Turning to the Higgs sector the vertex correction AC£^ is given by 

- (^D/j ^^^^^^^ 

+ "^x° (Cf),^ (^^i?)^/. Co (rnjo , 

+ "^x° "^x.- (^x^),, (^««)"/^ (^'^^)<^/. ^0 

+ "^X.-"^X° (^S)^^ (C:R)aJi (A^rfL)aJi ^0 ' "^x" ' "'I, 

Coo (mJ_,mJo,m|J-i 



(D.4) 



• ml-, ml 



+ MC^n {C:R),j,{NdL)^j, 



(D.5) 
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AC^^ may be obtained via the substitution L ^ R. Finally the corrections to the neutral 
Higgs vertex may be written 



,t^{c2 (3) m, [Sf) (GS^),, (G,,),^. Co ml , 



'Xa \ « 



+ 



1 2 
di' dj 

2 _2 2 



IJ 

"^X^^xl (^*dR)pM i^dL)aJ, Co {^ll^^ll^^X 



Xa Xt \ X ^/ab 



2 2 

m _ , m- 



+ 4 



+ 4 5: 



{^x°r) „^ (^*dR)f3M iNdL)ajj Coo (m^o , m^o , m|J - ^ 

Cnn fm^ _ , _ , m~ ^ 

V Xa ' Xt ' "V 8 



.50 

■'x- i?;^^ y^dRJbli K^dDalj 



{CdR)bIi {CdLlal 



(D.6) 



The correction to the right-handed vertex may be obtained in a similar manner to AC^^ . 
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